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Using  the  renormalization  group  all  the  couplings  of  the  Standard  Model 
and  its  minimal  supersymmetric  extension  are  run  to  two  loops,  taking  full  ac- 
count of  the  Yukawa  sector.  It  is  found  that  in  the  standard  model  the  gauge 
couplings  fail  to  unify,  whereas  extending  the  model  through  supersymmetry 
achieves  unification.  Bounds  are  placed  on  the  top  quark  mass  by  requiring 
equality  of  the  bottom  and  r  Yukawa  couplings  at  the  scale  of  unification.  In  its 
simplest  form,  the  supersymmetric  model  has  a  degenerate  superparticle  spec- 
trum, and  for  MsuS^'  =  1  TeV,  and  M5  =  4.6  GeV,  one  finds  139  <  A/<  <  194 
GeV,  which  remarkably  satisfy  the  p— parameter  bound.  The  corresponding 
bounds  on  the  Higgs  mass  are  found  to  be  44  <  Mfj  <  120  GeV.  The  model  is 
then  coupled  to  supergravity  yielding  a  richer  superparticle  spectrum  and  mak- 
ing it  more  accountable  to  experiment.  This  model  has  the  attractive  feature 
that  the  electro-weak  symmetry  is  radiatively  broken.  For  the  special  case  in 
which  global  supersymmetry  breaking  arises  solely  from  soft  gaugino  masses, 
Mt  is  expected  to  be  less  than  ~  130  GeV.  A  higher  upper  bound  is  predicted, 
if  Ml)  is  at  the  lower  end  of  its  experimental  uncertainty. 
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CHAPTER  1 
INTRODUCTION 

In  the  last  few  years,  it  has  become  apparent,  using  the  ever  increasing  ac- 
curacy in  the  measurement  of  the  strong  coupUng,  that  supersymmetry  (SUSY) 
affords  an  elegant  means  to  achieve  gauge  coupling  unification  [1]  at  scales  con- 
sistent with  grand  unified  theories  [2]  (GUTs).  Whereas  in  the  standard  model 
(SM)  the  three  gauge  couplings  unify  ''two  by  two"  forming  the  "GUT  trian- 
gle," in  the  simplest  minimal  supersymmetric  extension  of  the  standard  model 
(MSSM),  these  gauge  couplings  spectacularly  unify  at  a  point  (within  the  ex- 
perimental errors  in  their  values).  Given  that  the  scale  of  unification  in  these 
models  is  generally  above  the  lower  bound  set  by  proton  decay,  the  so-called 
SUSY-GUTs  have  gained  increasing  interest.  Constraints  coming  from  Yukawa 
couphng  unification  in  supersymmetric  SU{5)  and  50(10)  models  can  be  used 
to  yield  interesting  predictions  for  various  low  energy  parameters  including  the 
top  quark  mass    [3,4]. 

This  analysis  employs  the  renormalization  group  (RG)  [5]  to  extrapo- 
late the  parameters  of  the  standard  model  and  of  its  minimal  supersymmetric 
extension  to  unexplored  scales  [6].  With  the  ever  increasing  precision  of  ex- 
periment, the  inclusion  of  two  loop  effects  is  crucial.  The  complete  Yukawa 
sector  contribution  is  also  included.  In  the  standard  model  case,  thresholds 
effects  are  implemented  in  both  a  simple  and  a  more  involved  method.  The 
results  of  these  two  methods  are  discussed.  Numerical  methods  are  used  to 
evolve  the  parameters  to  different  scales  using  the  3  functions  foimd  in  the 
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literature  [7],  and  the  results  are  plotted  for  representative  values  of  the  Higgs 
boson  and  top  quark  masses.  The  running  of  the  quark  and  lepton  masses  and 
of  the  Cabibbo-Kobayashi-Maskawa  (CKM)  angles  is  generally  given  through 
the  running  of  the  Yukawa  matrices  (even  to  one  loop).  In  this  work,  the  quark 
masses  and  CKM  angles  are  evolved  by  diagonalizing  the  Yukawa  matrices  at 
every  step  in  the  Runge-Kutta  method  used  in  solving  the  j3  functions.  Often 
it  is  assumed  that  the  contribution  of  the  Yukawa  couplings  matrix  is  given 
essentially  by  the  top  quark  Yukawa  since  it  is  much  larger  than  the  others. 
Sometimes  a  better  approximation  is  made  by  keeping  only  the  diagonal  en- 
tries. The  present  numerical  technique  represents  a  small  improvement  over 
these  methods. 

Chapter  2  consists  of  a  basic  introduction  to  renormalization  and  renor- 
malization  scheme  dependence. 

Chapter  3  addresses  the  standard  model.  A  review  of  initial  data  extrac- 
tion from  experiments  is  presented.  Many  excellent  reviews  may  be  found  in 
the  literature,  e.g.,  Maxciano's  [8]  or  Peccei's  [9].  The  determination  of  the 
standard  model  gauge  couplings  is  discussed  as  is  the  initial  data  extraction  of 
the  Yukawas  and  the  CKM  angles.  The  extraction  of  the  quark  masses  from 
data  is  also  discussed.  This  is  a  complex  issue  well  known  to  be  marred  by 
the  nonperturbative  nature  of  QCD.  Hence,  in  the  low  energy  regime,  the  pure 
QCD  three-loop  contribution  is  included  in  the  analysis  of  the  running  of  the 
quark  masses.  Initial  data  for  lepton  masses  follow  the  quark  discussion.  Then 
the  extraction  of  and  constraints  on  the  physical  top  and  Higgs  masses  are 
considered.  The  scale  dependence  of  the  renormalized  scalar  \'acuum  expecta- 
tion value  is  addressed.  The  method  used  to  obtain  the  values  of  all  running 
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parameters  at  the  same  initial  scale  is  described  in  Appendix  E.  In  the  fol- 
lowing section,  threshold  effects  are  discussed.  Finally,  a  quantitative  analysis 
of  the  results  is  presented.  The  effects  of  using  one  loop  versus  two  loop  /3 
functions  are  contrasted  and  of  including  a  proper  versus  a  naive  treatment  of 
thresholds.  Plots  of  all  the  running  parameters  over  the  entire  range  of  mass 
scales  are  included  and  also  used  to  display  the  effects  discussed.  Furthermore 
some  tables  are  presented  with  actual  numerical  differences  associated  with 
these  effects. 

Chapter  4  addresses  the  minimal  supersymmetric  extension  of  the  stan- 
dard model.  A  brief  discussion  of  supersymmetric  models  is  presented.  The 
procedure  by  which  top  quark  mass  bounds  axe  determined  using  the  equality 
of  the  bottom  and  tau  masses  at  the  scale  of  unification  is  discussed  as  well 
as  the  one  light  Higgs  limit  employed.  The  initial  data  for  gauge  couplings 
and  quark  masses  are  presented.  This  is  followed  by  results  on  the  top  quark 
and  Higgs  boson  masses.  Finally,  a  discussion  on  how  to  improve  the  results 
is  presented. 

Chapter  5  implements  some  of  the  improvements  discussed  at  the  end  of  the 
last  chapter.  Namely,  the  supersymmetric  two  loop  /3  functions  are  included. 
Also,  soft  symmetry  breaJcing  terms  are  added.  These  lead  to  a  nondegenerate 
superparticle  spectrum  as  well  as  to  the  radiative  breaking  of  the  electro- weak 
symmetry.  Similar  analyses  have  appeared  in  the  literature  but  use  one  loop 
/?  functions  and  the  tree  level  Higgs  potential  [10].  In  Chapter  5,  a  brief 
discussion  of  the  effective  one  loop  potential  is  presented.  The  mass  formulas 
for  the  sfermions.  Higgses,  charginos,  and  neutralinos  are  given.  The  boundary 
conditions  at  the  unification  scale  in  these  minimal  low  energy  supergravity 
models  are  discussed.    Ne.xt  the  numerical  procedure  employed  is  described. 


4 
The  treatment  of  thresholds  and  the  "special"  form  of  the  l3  functions  needed 
is  discussed  next.  Finally,  some  preliminary  results  are  presented. 

Chapter  6  contains  the  conclusions  of  this  work. 

Appendix  A  contains  all  the  needed  two  loop  l3  functions  for  renormaliza- 
tion  group  studies  of  the  standard  model.  Appendix  C  contains  those  of  the 
minimal  supersymmetric  extension  of  the  standard  model.  Many  of  these  /3 
functions  have  yet  to  appear  in  the  literature  in  as  general  a  form.  Appendix 
B  presents  some  examples  in  calculation  of  the  /?  functions  of  Appendix  C. 
Appendix  D  deals  with  issue  of  the  vacuum  expectation  value's  l3  function.  A 
toy  model  is  used  to  gather  insight  into  the  problem.  The  numerical  solution 
routines  used  extensively  throughout  this  work  are  discussed  in  Appendix  E. 
Appendix  F  contains  a  cumbersome  formula  needed  in  the  extraction  of  the 
Higgs  boson  physical  mass.  Finally  Appendix  G  presents  a  novel  use  of  the 
spinor  representations  of  the  orthosymplectic  Lie  supergroup  0Sp{r/2m,R) 
[11]. 


CHAPTER  2 
THE  RENORMALIZATION  GROUP 

Renormalization  is  a  reparametrization  of  a  theory  which  renders  Green 
functions  and  physical  quantities  finite  order  by  order  in  perturbation  theory. 
A  specific  choice  of  renormahzed  parameters  defines  a  renormahzation  scheme. 
The  physics  is,  of  course,  independent  of  how  the  theory  is  renormahzed.  A 
common  way  of  relating  bare  and  renormahzed  parameters  is 

9o  =  9-Sg,  (2.1) 

where  go  is  the  bare  parameter,  g  is  the  renormahzed  parameter,  and  6g  is 
the  counterterm.  Fixing  the  counterterms  by  requiring  them  to  consist  only 
of  the  infinite  terms  needed  to  render  the  theory  finite  defines  the  minimal 
subtraction  (MS)  prescription  [12].  A  feature  of  the  MS  scheme  is  a  mass 
scale  ^  which  enters  in  the  process  of  regularizing  divergent  integrals  using 
dimensional  regularization.  Furthermore  the  unit  of  mass  fi  is  used  to  keep 
couplings  dimensionless  when  continuing  to  d  dimensions  in  the  dimensional 
regularization  procedure.  For  example,  if  Eq.  (2.1)  represents  any  of  the  three 
gauge  couplings  of  the  Standard  Model,  then  /u  is  introduced  as  follows  to  keep 
them  dimensionless 

9o{Qti)~' =  g  -  Sg  ,  (2.2) 

where  ^  is  a  constant  parametrizing  the  arbitrariness  in  the  finite  parts  of  diver- 
gent integrals  in  dimensional  regularization  and  e  =  {4  —  d)/2.  Equation  (2.2) 
defines  a  family  of  MS  schemes.    Choosing  ^  =  1  is  the  simplest  MS  scheme 
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which  was  described  above.    Choosing  ^  =  e"^^ /An,  where  7^  =  0.5722... 
is  the  Euler-Mascheroni  constant,  defines  the  so-called  modified  minimal  sub- 


traction (MS)  prescription  [13].  This  scheme  is  the  most  commonly  employed 
in  QCD  calculations,  and  it  is  the  one  adopted  here.  The  free  parameters  of 
the  Standard  Model  in  the  MS  schemes  Eire  /i  dependent.  Their  /i  evolution 
is  governed  by  the  0  functions  of  the  renormalization  group.  Moreover,  these 
running  paraxneters  are  not  in  general  equal  to  their  corresponding  physical 
values  (consequently,  for  the  masses,  a  convention  is  adopted  wherein  upper 


case  M  refers  to  physical  values  and  lower  case  m  denotes  MS  values).  This 
is  to  be  contrasted  with  the  on-shell  renormzJization  scheme  in  which,  for 
example,  the  renormalized  masses  equal  their  physical  values  and  the  renor- 
malized  electromagnetic  coupling  equals  the  fine  structure  constant.  However, 
the  MS  schemes  have  the  attractive  characteristic  that  the  f3  functions  are  /i 
independent  and  therefore  particularly  simple  to  integrate.  Physical  quanti- 
ties P{{gi{fj,)},fi)  expressed  in  terms  of  /x  and  the  running  parameters  of  the 
theory,  {giifJ.)},  must  be  /z  independent 

fi  —  P{{giip)},^i)  =  (^—  +  0i  —  )P  =  0  ,  (2.3) 

where  the  0i  are  the  3  functions.  The  two-loop  l3  functions  of  the  Standard 
Model  have  been  collected  in  Appendix  A. 

As  mentioned  above  physical  quantities  are  renormalization  scheme  inde- 
pendent. However,  this  assumes  that  calculations  can  be  done  without  ap- 
proximation. In  reality,  calculations  are  only  perturbative  approximations  and 
these  do  depend  on  the  renormalization  scheme. 

Consider  massless  QCD  with  the  one  dimensionless  coupling,  as-  Suppose 
a  physical  quantity,  P.  is  calculated  to  n      order  in  perturbation  theory  in  two 
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renormalization  schemes,  then  the  n      order  approximation  is  given  by 


Pn  =  asipo  +  Piois  +  P20t'i -\ +  Pn-l^l    ^)  ,  (2.4) 


in  one  scheme  and  by 

Pn=0's{Po+PlO's+P2'^s     +---+Pn-l«s         ),  (2.5) 

in  the  other  scheme.  The  two  couphngs,  defined  in  their  respective  renormal- 
ization schemes,  can  be  related  to  each  other 

a',  =  as{l  +  bias  +  b2al  +  ■■■  +  ftn-ia^^)  .  (2.6) 

Substituting  Eq.  (2.6)  into  Eq.  (2.5)  yields 

d'  I   II    t      II         ,      II    2    ,  ,11  n—l\ 

Pn=(^s(Po+ PIG'S +P20is  +  ---+Pn-l(^s       I 

(2.7) 
+  a,(pM  +  ---+P2n-l«?""')- 

As  this  is  an  approximation  to  the  same  quantity,  P,  the  first  n  terms  can  be 

identified  with  Pn-  Therefore  the  two  approximations  differ  by  terms  of  higher 

order 

P^  -  Pn  =  a,(p'^a^  +  •  •  •  +  p^'^.ia?"-!)  .  (2.8) 

In  QCD  where  Og  is  large,  this  difference  may  be  large  and  thereby  lead  to 
renormalization  scheme  dependence  problems.  In  the  Standard  Model  as  in 
QED  where  the  couplings  are  small,  this  is  not  so  great  a  problem.  In  QCD 
where  the  strong  couphng  as  is  large,  there  will  be  renormahzation  scheme 
dependence  problems.  In  the  electroweak  model  as  in  QED  where  the  couplings 
are  small,  this  is  not  so  great  a  problem. 


CHAPTER  3 
THE  STANDARD  MODEL 


3.1  a\iMz)^  a'^jMz).  and  a-}(Mz) 
The  determination  of  the  SU(2)l  x  U(1)y  couphngs  proceeds  from  the 
Standard  Model  relations 


ai(/^)  =  -7 —  =  C  2d     t    \ 

47r  cos''  (^w{fi) 


2,     ^  .      s  (3.1.1) 


a  (    )=  ^2^^^  _        Q'(/^) 


47r  sin2^py(/x) 

where  af/x)  =  e  (/x)/47r  and  C  is  a  normalization  constant  which  equals  1 
for  the  Standard  Model  and  equals  |  when  the  Standard  Model  is  incorpo- 
rated in  grand  unified  theories  of  the  SU{N)  and  SO{N)  type  [5].  What 
is  required  to  specify  these  couplings  are  the  values  of  a(^)  and  sin  d\^(fi) 
in  the  renormalization  scheme  employed  (i.e.,  MS).  The  electromagnetic  fine 
structure  constant  (a^m  ~  137.036)  is  extrapolated  from  zero  momentum  scale 
to  a  scale  fj,  equal  to  AI^  in  the  present  case.  In  pure  QED  with  one  species 
of  fermion  with  mass  m  the  MS  renormalized  vacuum  polarization  function  is 
given  by 

2  ^  2 

nC?^)  =  4^(ln^  -  6  / dx  x(l  -  x)ln[l  -  x(l  -  x)\])  .  (3.1.2) 

OK         m'^         J  m^ 

0 

The  renormalized  coupling  a(iJ.)  is  related  to  the  fine  structure  constant  aem 

as  follows 

a(/i) 


OCpm   — 


i  +  n(0) 
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(3.1.3) 
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In  the  Standard  Model  where  there  are  many  species  of  charged  fermions  and 
charged  gauge  bosons,  Eq.  (3.1.3)  generahzes  to    [14] 

^-H/^)  =  «e"m  -  :r  E  ^/^(^  -  m;)ln-^  +  ^  .  (3.1.4) 

OTT  '—-'       •'  '  nif  DTT 

The  effects  of  the  strong  interaction,  which  enter  as  a  hadronic  contribution  to 
the  vacuum  polarization  function,  must  be  included  also.  The  nonperturbative 
nature  of  the  strong  interaction  at  low  momentum  is  handled  by  rewriting  the 
hadronic  contribution  to  the  vacuum  polarization  at  zero  momentum  as 

n''(0)  =  (n'^(o)  -  "^^{s^))  +  n'*(g2) .  (3.1.5) 

If  q  is  chosen  large  enough,  11  (g)  can  be  calculated  perturbatively.  The 
terms  (H  (0)  — 11  (5  ))  can  then  be  related  to  the  total  cross  section  for  e^ e~  — + 
hadrons    [14].  Using  the  optical  theorem,  one  may  write 

Im{n''(5)}  =  — ^ — cr(e+e~  -^  hadrons)  ,  (3.1.6) 

47raem 

where  s  is  the  square  of  the  center  of  mass  energy.  For  the  process  e^e~  -^ 
\i^ pT ,  the  cross  section  is  calculated  to  be  (taking  va^  =  0) 

a(e+e--./z+^-)  =  ^^.  (3.1.7) 

3s 

In  terms  of  the  ratio  of  these  two  cross  sections. 


a{e'^e     —>■  hadrons) 
a-(e+e~  — >  n'^fi' 


Ris)  =  "^    ;    -—-7  ,  (3.1.8) 


one  may  write  Eq.  (3.1.6) 


Imjn^s)}  =  ^R{s)  .  (3.1.9) 
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Using  an  unsubtracted  dispersion  relation  for  11  (q)  the  combination  (11  (0) — 
U.  {q  ))  can  be  expressed  as 

n'm  -  ^\f)  =  ^  /  rf.  -1^  .  (3.1.10) 

4ml 
This  can  be  evaluated  using  experimentally  known  data.  This  procedure  yields 
a  value 

a-\Mz)  =  127.9  ±  0.3  .  (3.1.11) 

The  process  independent,  renormalized  weak  mixing  Eingle  sin   d\Y  of  the 
on-shell  scheme  is  defined  to  be 

sin2%  =  l ^,  (3.1.12) 

M| 

where  My[r  and  Mz  are  the  physical  masses  of  the  W  and  Z  gauge  bosons. 
Knowing  the  precise  values  of  the  W  and  Z  boson  masses  and  using  the  equa- 
tion above  provide  one  way  of  extracting  the  value  of  sin  9\y.  Alternatively, 
the  bare  relation  involving  the  low  energy  Fermi  constant  measured  in  muon 
decay  and  the  W  boson  mass 


G,o  e± 


\/2        Ssin2%„Mpy2  ' 
may  be  corrected  to  order  a  and  rewritten    [15.16] 


(3.1.13) 


Mw  =  Mz  cos%  =  (^^)' ^- r  ,  (3.1.14) 

v-^Crp      sin^py(l  -  Ar)2 

with  [T^oteml v^Gpi)^  =  37.281  GeV  and  Ar  is  a  parameter  containing  order 
a  radiative  corrections  and  which  depends  on  the  mass  of  the  top  and  Higgs. 
The  radiative  corrections  represented  by  Ar  can  be  viewed  as  accounting  for 
the  mismatch  in  the  scales  associated  with  the  parameters  of  the  relation.  G^ 
and  atm  are  low  energy  parameters  whereas  M\^  and  sin   6\y  axe  associated 
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with  the  electroweak  scale.  One  can  absorb  the  radiative  effects  using  the 
renormahzation  group  by  replacing  G^  and  aem  with  corresponding  running 
parameters  at  M^ 

""^^^)  ~1.  (3.1.15) 


v/2G^(M^)Af2,sin2% 
Combining  Eqs.  (3.1.14)  and  (3.1.15)  gives 

aem     GuiMz)  . 

Arsil TTJ— 7— T^; •  (3.1.16) 

a{Mz)      G/i 

Using  Eq.  (3.1.4)  and  the  fact  that  G^{Mz)  ~  G^  (see  Section  2.8)  gives  an 
estimate  of  the  size  of  the  radiative  corrections 

Ar«0.07.  (3.1.17) 

For  large  values  of  Mt  and  Mfj  {Mt,  Mfj  >  Mz)    [15,17] 

A.  «  1  -  -^  -        ^"7       4  +       ''"T      In^  .  (3.1.18) 

a(Mz)       167rsin''%M|       487rsin2%      M| 

A  third  way  of  extracting  sin  B^r  is  from  neutral  current  experiments,  among 
which  deep  inelastic  neutrino  scattering  appears  to  provide  the  best  determi- 
nation. A  running  sin  ^^(/x)  may  be  defined  in  MS  and  differs  from  the  above 
sin^  ^j,y  by  order  a  corrections.  The  MS  running  W  boson  mass  m^,^(/i)  and 
the  corresponding  physical  mass  M\.y,  identified  as  the  simple  pole  at  q  —  My^ 
of  the  W  propagator,  are  related  as  follows 

M'lr  =  n?^(^i)  +  Aj^^y(M^,^)  ,  (3.1.19) 

where  Ai^yy^  is  the  transverse  part  of  the  W  self-energy.  A  similar  relation 
holds  for  the  Z  boson.  In  MS  renormahzation,  the  following  relation  defines 
the  running  sin   9\y{^) 

sin2%(/x)  =  l--^  .  (3.1.20) 

m|(//) 
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Equation  (physmw)  and  its  Z  analog  may  be  combined  with  Eq.(3.1.20)  to 

give 

sin^%(/^)  ^  ^  _  cos2%^4^(M|,/x)  _  Af^^jM^,,^)^  ^^^^^^^^ 

sin^^py  sin^^py  M|  M^ 

An  explicit  expression  relating  sin^  By/  and  sin^  ^py(A/iy)  is  given  in  Ref.  18. 

Another  relation  for  sin^  ^w(^)  ^^7  ^^  arrived  at  directly  linking  it  to  M^ 

[19]  or  Myi/     [20].    In  paxticulax,  if  one  chooses  A/py  as  the  input  mass,  then 

one  introduces  a  radiative  correction  parameter  Arpy  such  that 

sin^  %(M^)(1  -  Afpy)  =  sin^  %(1  -  Ar)  ,  (3.1.22) 

from  which  it  follows  that 

sin2%(M2)  =      I     ■"     '         .  (3.1.23) 

M^(l  -  Arpy) 

Similarly  one  can  introduce  a  radiative  correction  Ar^^  if  one  chooses  M^  as 
the  input  mass 

sin^  0w{Mz)  cos'^  %(M^)(1  -  Afz)  =  sin^  %  cos^  9^(1  -  Ar)  .     (3.1.24) 

A  fit  to  all  neutral  current  data  gives 

sin^  OwiMz)  =  0.2324  ±  0.0011  ,  (3.1.25) 

for  arbitrary  Mf    [21].  Using  these  values  of  a{AIz)  and  siiir  6^r(Mz)  yields 

ai{Mz)  =  0.01698  ±  0.00009  , 

(3.1.26) 
a2{Mz)  =  0.03364  ±  0.0002  . 

The  value  of  the  strong  coupling  has  not  been  determined  experimentally  as 
well  as  ai  and  02.  The  nonperturbative  nature  of  low  energy  QCD  leads  to  rel- 
atively large  uncertainties  in  its  determination.  Different  processes  for  the  ex- 
traction of  0:3  yield  significantly  different  results.  Unless  otherwise  stated,  the 
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value  of  03  used  in  the  runs  will  be  that  obtained  from  a  Gaussian  weighted  av- 
erage of  the  results  of  several  processes,  including  e"'"e~  scattering  into  hadrons 
[22],  heavy  quarkonium  decay    [23],  deep  inelastic  scattering    [24],  and  e+e" 
scattering  into  jets    [25].  The  value  found  is    [6] 

a3(M2)  =  0.113  ±0.004  .  (3.1.27) 

3.2  Yukawas 

To  take  full  account  of  the  Yukawa  sector  in  running  all  the  couplings,  initial 
\Tilues  for  the  Yukawa  couplings  are  necessary.  They  must  be  extracted  from 
physical  data  such  as  quark  masses  and  CKM  mixing  angles.  Furthermore,  the 
interesting  parameters  to  be  plotted  must  be  determined  step  by  step  in  the 
process  of  running  to  Planck  mass.  These  two  procedures  are  not  unrelated 
and  require  the  diagonalization  of  the  up-type,  down-type,  and  leptonic  Yukawa 
matrices. 

Machacek  and  Vaughn's  [7]  convention  are  used  where  the  interaction 
Lagrangian  for  the  Yukawa  sector  is 

C  =  uYu^^Q  +  dYj^^Q  +  eYe^^L  +  h.c.  .  (3.2.1) 

The  Yukawa  couplings  are  given  in  terms  of  3x3  complex  matrices.  After  elec- 

troweak  symmetry  breaking,  these  translate  into  the  quark  and  lepton  masses 

^9  /me        0         0    \ 
Ye  =  —        0        m^        0  , 

^'     \    0  0        mrj 

s/2  f^d       0        0   \ 
Yd  =  —        0       ms       0         ,  (3.2.2) 


^     ^    0        0       m 


h. 


J^  (mu      0        0  \ 
Y«  =  —       0       mc      0       y  , 


^'       V      0  0  777 
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where  V  is  the  CKM  matrix  which  appears  in  the  charged  current 

J^  ~  ULlfiVdL  .  (3.2.3) 

It  is  a  unitary  3x3  matrix  often  parametrized  as  follows 

(ci  51C3  -S1S3  \ 

-51C2     C1C2C3  -  52536"^     C1C253  +  52036*^        ,  (3.2.4) 

-51^2     ciS2C^  +  C2S^e^^     cis^s^  -  C2Czt^^  j 

where  si  =  sin^,-  and  Ci  =  cos^j,  i  =  1,2,3. 

The  entries  of  the  parametrized  CKM  matrix  can  be  related  simply  to 

the  experimentally  known  CKM  entries.    The  paxticle  data  book     [21]  gives 

the  following  ranges  of  values  (assuming  unitarity)  for  the  magnitudes  of  the 

elements  of  the  CKM  matrix 

/ 0.9747-0.9759       0.218-0.224         0.001-0.007   \ 
\V\  =        0.218-0.224       0.9734-0.9752       0.030-0.058         .  (3.2.5) 

V   0.003-0.019        0.029-0.058      0.9983-0.9996/ 

These  ranges  of  values  can  be  converted  to  bounds  for  Si,  i  =  1, 2, 3,  and  sin^. 

These  bounds  are  arrived  at  by  finding  values  for  the  four  angles  such  that  the 

entries  of  the  CKM  matrix  obtained  from  these  satisfy  the  conditions  imposed 

by  Eq.(3.2.5).  One  finds 

0.2188  <  sin  ^1  <  0.2235  , 

0.0216  <  sin  ^2  <  0.0543  ,  (3.2.6) 

0.0045  <  sin  ^3  <  0.0290  . 
However,  the  accuracy  with  which  \V\  is  known  does  not  constrain  sin^.  A  set 
of  angles  {^i,^25%?<^}  was  chosen  that  falls  within  the  ranges  quoted  above. 
The  initial  data  needed  to  run  the  Yukawa  elements  are  extracted  from  the 
CKM  matrix  and  the  quark  masses.  A  problem  arises  though  for  the  mixing 
angles,  which  was  solved  for  the  quark  masses  (see  Section  2.5),  in  that  it 
is  not  clear  at  what  scale  the  chosen  initial  values  for  these  angles  should 
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be  considered  known.  However,  it  is  observed  that  for  the  whole  range  of 
initial  values,  the  running  of  the  mixing  angles  is  quite  flat,  with  a  perceptible 
increase  in  62  between  M\v  and  the  Planck  scale  for  higher  top  masses.  This 
is  in  accordance  with  the  angles  being  related  to  ratios  of  quark  masses,  and 
therefore,  the  exact  knowledge  of  that  scale  (or  lack  thereof)  is  not  as  critical 
as  might  be  feared  a  'priori. 

3.3  Quark  and  Lepton  Masses 

There  are  large  theoretical  uncertainties  in  the  extraction  of  the  masses 

of  the  three  Hghtest  quarks  from  experiment.   They  are  determined  by  chiral 

perturbation  techniques  and  QCD  spectral  sum  rules   [26.27].  In  the  following, 

the  MS  running  masses  at  1  GeV  of  the  three  lightest  quark  will  be  taken  to 

be    [6] 

mu{l  GeV)  =  5.2  ±  0.5  MeV  , 

m^(l  GeV)  =  9.2±0.5  MeV  ,  (3.3.1) 

ms(l  GeV)  =  194  ±4  MeV  , 
For  the  charm  and  bottom,  the  nonrelativistic  bound  state  approximation  may 

be  applied.  One  speaks  of  physical  masses  and  associates  these  with  the  pole  of 

the  quark  propagator.  A  weighted  average,  based  on  results  from  J/t/'  and  T 

sum  rules    [28]  and  from  heavy-light,  B  and  B* ,  D  and  D*  meson  masses  and 

semileptonic  B  and  D  decays    [29],  yields  the  following  physical  and  running 

masses  at  1  GeV  for  the  charm  and  bottom    [6] 

Mc  =  1.60  ±  0.05  GeV,     mc(l  GeV)  =  1.41  ±  0.06  GeV  , 

(3.3.2) 
M^  =  4.89  ±  0.04  GeV,     m,,{l  GeV)  =  6.33  ±  0.06  GeV  , 

Although  the  associated  errors  on  these  averages  are  restrictive,  there  is  larger 

uncertainty  in  the  actual  central  values.    For  example,  in  the  bottom  quark. 

values  as  small  as  ?77j(l  GeV)  =  5.7  GeV  and  as  large  as  mft(l  GeV)  =  6.5 
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GeV  are  acceptable.  In  the  case  of  the  Ughter  quarks,  their  mass  ratios  are 
known  more  accurately  than  the  actual  value  of  their  respective  masses.  The 
ratios  m^lmu  =  1.8  and  mslm^  =  21  are  widely  accepted. 

The  physical  (pole)  masses  of  the  leptons  are  very  well  known    [21] 

Me  =  0.51099906  ±  0.00000015  MeV  , 

Mp  =  105.658387  ±  0.000034  MeV  ,  (3.3.3) 

Mr  =  1.784llg:gg27  GeV  . 

These  masses  represent  an  example  of  a  physical  quantity  as  discussed  in  Chap- 
ter 2.  Indeed  expressed  in  terms  of  the  renormalized  parameters  (i.e.,  masses 
and  couplings)  of  the  theory,  the  physical  meiss  is  just  the  simple  pole  of  the  rel- 
evant field's  propagator.  Suppose  S{^;g,m,ij.)  is  the  renormalized  propagator 
of  some  fermion  with  a  simple  pole  at  ^  =  M{g,  m.,  fi)  such  that 

S(ff;  g,  m,  n)  =  - — — ,  (3.3.4) 

p  -  M(g,m,ii) 

with  finite  residue,  z{j^  =  M;g,m,^).  The  relevant  renormalization  group 
equation  is 

ifi—  +  0-Q-  +  "^^rn-g^  +  1-f)S{j};  g,m,fi)  =  0  ,  (3.3.5) 

where  7  is  the  anomalous  dimension  and  jm  is  the  mass  anomalous  dimension. 
Inserting  Eq.  (3.3.4)  in  Eq.( 3.3.5),  then  multiplying  the  resulting  expression  by 
(^  —  M)  ,  and  lastly  setting  jf  —  M  gives  the  renormalization  group  equation 
for  M 

H—M{g,m,n)  =  {fi—  +  ^—  +  mfm-g^)M{g,m.,  fi)  =  0  .  (3.3.6) 

These  values  are  used  to  determine  initial  data  for  the  running  masses.  Some 
authors  neglect  QED  corrections  and  use  the  physical  values  for  the  running 
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values  at  ~  Mz-,  which  introduces  only  a  small  error.  By  calculating  the  one- 
loop  self-energy  corrections,  one  arrives  at  a  QED  relation  between  the  running 
MS  masses  and  the  corresponding  physical  masses 

mM  =  Mi[\  -  ^(In^  +  \)]  .  (3.3.7) 

Choosing  /i  =  1  GeV  as  in  the  quark  mass  case  and  using  Eqs.  (3.3.7)  and 
(3.1.4)  yields  the  running  lepton  masses  (taking  m/  =  M/  in  the  log  term 
above  is  an  appropriate  approximation  to  order  a) 

me(l  GeV)  =  0.4960  MeV  , 

m^(l  GeV)  =  104.57  MeV  ,  (3.3.8) 

mr(l  GeV)  =  1.7835  GeV  . 

3.4  Top  and  Higgs  Masses 
The  Higgs  boson  and  top  quark  masses  have  not  been  measured  directly  at 
present;  however,  their  values  affect  radiative  corrections  such  as  Ar.  Consis- 
tency with  experimental  data  on  sin  9]/i/  requires  Mt  <  197  GeV  for  M//  =  1 
TeV  at  99%  CL  assuming  no  physics  beyond  the  Standard  Model  [30].  Pre- 
cision measurements  of  the  Z  mass  and  its  decay  properties  combined  with 
low  energy  neutral  current  data  have  been  used  to  set  stringent  bounds  on 
the  top  quark  mass  within  the  minimal  Standard  Model.  A  global  analysis  of 
this  data  yields  Mt  =  I22II2  GeV,  for  all  allowed  values  of  jU//  [31].  Recent 
direct  search  results  set  the  experimental  lower  bound  Mt  ^  91  GeV.  As  for 
the  Higgs,  the  analysis  of  Ref.  31  gives  the  restrictive  bound,  Mf{  ^  600  GeV, 
if  Mt  <  120  GeV,  and  Mf{  <  6  TeV.  for  all  allowed  ^/^  Since  perturbation 
theory  breaks  down  for.  A///  ^  1  TeV.  the  latter  bound  on  the  Higgs  boson 
mass  is  not  necessarily  meaningful.  LEP  data  set  a  lower  bound  on  the  Higgs 
boson  mass  of  48  GeV    [32]. 
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In  the  present  analysis,  initial  values  of  the  MS  running  top  quark  mass 
mt  and  of  the  scalar  quartic  self-coupling  A  at  M^  are  chosen  arbitrarily  (con- 
sistent with  the  bounds  quoted  above).  As  noted  earlier  in  Chapter  2,  these 
running  parameters  are  not  equal  to  their  physical  counterparts.  However,  any 
reasonable  prediction  for  the  masses  of  the  top  quark  and  of  the  Higgs  boson 
that  may  come  from  this  analysis  should  be  that  of  experimentally  relevant, 
physical  masses.  Therefore,  formulas  similar  to  Eq.  (3.3.7)  relating  MS  running 
parameters  to  physical  masses  axe  needed.    To  calculate  the  physical  or  pole 

mass  of  the  top  quark,  the  following  equation  is  used    [33] 

5 

—-{—  =  1  +  -      '     "  +  16.11  -  1.04  y(l  -  -rj-  )    ,      (3-4.1 

mt{Mt)  3       K  ^'        M<  ''       TT       '    '      ' 

where  A/f,  i  =  1, . .  .  ,  5,  represent  the  masses  of  the  five  lighter  quarks.  Like- 
wise the  physical  mass  of  the  Higgs  boson  can  be  extracted  from  the  following 
relation    [34] 

A(//)  =  ^mJ(1  +  6(^)),  (3.4.2) 

where  S{fj,)  contains  the  radiative  corrections.  Its  form  is  rather  elaborate  and 
it  is  relegated  to  Appendix  B.  Equations  (3.4.1)  and  (3.4.2)  are  highly  nonlin- 
ear functions  of  Alt  and  Mff,  respectively.  Their  solution  requires  numerical 
routines  that  are  described  in  Appendix  C. 

At  the  one  loop  level,  the  gauge  couplings  are  unaffected  by  the  other 
couplings  in  the  theory.  On  the  other  hand,  the  Yukawa  couplings  are  affected 
at  one  loop  by  both  the  gauge  and  Yukawa  couplings.  Since  the  top  Yukawa 
coupling  is  at  least  as  big  as  the  gauge  couplings  at  low  energy,  that  means 
the  running  of  the  Yukawas  is  sensitive  to  mostly  the  top  Yukawa  and  the 
QCD  gauge  couplings.  Thus  one  expects  the  mass  and  mixing  relations  just 
described  to  be  sensitive  to  the  value  of  the  top  quark  mass.  The  Higgs  quartic 
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self-coupling  enters  in  the  running  of  the  other  couplings  only  at  the  two  loop 
level,  so  that  its  effect  on  the  quark  and  lepton  parameters  is  small.  However, 
its  own  running  is  very  sensitive  to  the  top  quark  mass;  it  can  become  negative 
as  easily  as  it  can  blow  up,  corresponding  to  vacuum  instability  or  to  strong  self- 
interaction  of  the  Higgs  (triviality  bound),  respectively  [35].  The  discovery 
of  the  Higgs  with  mass  outside  these  bounds  would  be  a  signal  for  physics 
beyond  the  Standard  Model.  The  graphs  in  Figs.  1-4  summarize  these  bounds 
for  representative  values  of  the  top  quark  mass.  For  example,  if  Mt  =  150  GeV, 
one  can  see  from  the  corresponding  plot  that  a  Higgs  mass  between  95  and  150 
GeV  need  not  imply  any  new  physics  up  to  Planck  scale.  However,  if  the  Higgs 
were  observed  outside  of  this  range,  then  some  new  physics  must  appear  at  the 
scale  indicated  by  the  curve,  either  because  of  vacuum  instability  if  Mfj  <  95 
GeV  or  because  the  Higgs  interaction  becomes  too  strong  if  Mjj  >  150  GeV.  It 
is  amusing  to  note  that  it  is  for  comparable  values  of  the  top  and  Higgs  masses 
that  these  bounds  are  least  restrictive,  but  it  is  important  to  emphasize  that 
a  high  value  of  the  top  with  a  relatively  low  value  of  the  Higgs  necessarily 
indicates  the  presence  of  new  physics  within  reach  of  the  SSC. 

3.5  Vacuum  Expectation  Value 
The  vacuum  expectation  value  (vev)  of  the  scalar  field  may  be  extracted 
from  the  well  known  lowest  order  relation 

t;  =  (\/2Gp)"  2  =246.22  GeV  .  (3.5.1) 

From  the  very  well  measured  value  of  the  muon  lifetime,  r^  =  2.197035  ± 
.000040  X  10~"  s  [21],  the  Fermi  constant  can  be  extracted  using  the  following 
formula   [36] 


where 


giving 
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fix)  =  l-8x  +  8x^  -x^  -  12x^lnx  ,  (3.5.3) 


Gf,  =  1.16637  ±  0.00002  x  10"^  GeV-^  .  (3.5.4) 


This  parameter  may  be  viewed  as  the  coefficient  of  the  effective  four-fermion 
operator  for  muon  decay  in  an  effective  low  energy  theory 

^[I7e7/?(1  -  75)e][/l7/?(l  "  IbM  ■  (3.5.5) 

A  direct  calculation  (e.g.,  in  the  Landau  gauge)  of  the  electromagnetic  correc- 
tions yields  that  the  operator  is  finitely  renormalized  (i.e.,  G^  does  not  run) 
[16,37].  Another  way  to  see  this  is  by  using  a  Fierz  transformation  to  rewrite 
the  above  expression 

^[i7e7^(l  -  I5>fi][lllp{l  -  75)e]  ■  (3.5.6) 

The  neutrino  current  does  not  couple  to  the  photon  field,  and  the  e  —  ^  current 
is  conserved  and  is  hence  not  multiplicatively  renormalized. 

An  initial  value  is  needed  for  the  running  vacuum  expectation  value  at  some 
scale  jj,.  Wheater  and  Llewellyn  Smith  [38]  consider  muon  decay  to  order  a 
in  the  context  of  the  full  electroweak  theory  and  derive  an  equation  relating 


an  MS  running  G^  to  the  experimentally  measured  value.  From  this  formula 
one  can  extract  a  value  for  v(Mz)-  However,  the  formula  is  derived  in  the  't 
Hooft-Feynman  gauge,  and  the  evolution  equation,  Eq.  (L18)  of  Appendix  A 
for  the  vev,  is  valid  only  in  the  Landau  gauge  (see  Appendix  D).  Nevertheless, 
motivated  by  the  discussion  of  the  previous  paragraph,  the  initial  condition  for 
the  vev  is  chosen  to  be  v{Miy)  =  246.22  GeV.  Using  the  initialization  algorithm 
(see  Appendix  C),  one  arrives  at  u(A/^).  It  is  found  that  this  procedure  leads 
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to  no  significant  correction,  and  one  therefore  takes,  ab  initio,  v{Mz)  =  246.22 
GeV. 

3.6  Thresholds 

In  mass  independent  renormalization  schemes,  the  running  coupHngs  are 
unphysical.  From  the  decoupHng  theorem  [39]  one  expects  the  physics  at  en- 
ergies below  a  given  mass  scale  to  be  independent  of  the  particles  with  masses 
higher  than  this  threshold.  Therefore,  for  a  correct  interpretation  of  these  run- 
ning couplings,  one  must  take  into  account  the  thresholds  [40,41.42].  For  the 
electroweak  threshold,  one  loop  matching  functions  [42]  are  used  with  the  two 
loop  beta  functions  valid  in  the  Standard  Model  regime  below  the  SUSY  scale. 
These  matching  functions  are  obtained  in  MS  renormalization  by  integrating 
out  the  heavy  gauge  fields  in  such  a  way  that  the  remaining  effective  action  is 
invariant  under  the  residual  gauge  group  [41].  At  the  electroweaJc  threshold, 
near  M\y,  the  heavy  gauge  fields  and  the  top  quark  axe  integrated  out.  Below 
this  threshold  there  is  an  effective  SU{3)c  x  U{l)£\f  theory.  Thresholds  in 
this  region  are  obtained  by  integrating  out  each  quark  to  one  loop  at  a  scale 
equal  to  its  physical  mass.  At  these  scales  the  one  loop  matching  functions 
in  the  gauge  couplings  vanish  and  the  threshold  dependence  appears  through 
steps  in  the  number  of  quark  flavors  [43]  as  the  renormalization  group  scale 
passes  each  physical  quark  mass. 

In  the  SM  runs  (see  next  section),  in  which  an  analysis  of  the  relative 
importance  of  including  a  proper  treatment  of  thresholds  effects  instead  of 
using  a  simple  step  function  technique,  it  is  found  that  the  former  method 
does  not  improve  significantly  over  the  latter   [6]. 
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3.7  Analysis  and  Results 
The  results  of  numerically  integrating  the  3  functions  of  the  Standard 
Model  parameters  from  1  GeV  to  Planck  mass  are  depicted  in  the  following 
figures.  For  most  of  these  plots,  the  arbitrary  choice,  Mt  =  Mf{  =  100  GeV,  is 
made.  Figure  5  displays  the  evolution  of  the  inverse  of  each  of  the  three  gauge 
couplings  including  the  associated  uncertainties  in  their  values.  In  it,  one  sees 
the  ''GUT  triangle"  signifying  the  absence  of  grand  unification,  assuming  the 
Standard  Model  as  an  effective  theory  in  the  desert  up  to  the  Planck  scale. 
Here,  the  differences  between  one-  and  two-loop  evolution  appear  in  the  high 
energy  regime.  Differences  are  also  present  for  the  strong  coupling  at  low  en- 
ergies where  it  becomes  large.  Note  that  the  uncertainties  do  not  fill  in  the 
"GUT  triangle."  Figures  6,  7,  and  8  display  the  evolution  of  the  light  mass 
fermions  (me,  m^,  and  rrid),  the  intermediate  mass  fermions  (m^  and  mg),  and 
the  heavy  mass  fermions  (rrz,-,  rric,  and  rrif,),  respectively.  The  largest  differ- 
ences between  one-loop  and  two-loop  evolution  occur  in  the  bottom,  charm, 
and  strange  quark  masses  in  these  cases.  In  Fig.  9,  the  quartic  self-coupling  A 
and  the  top  Yukawa  coupling  yt  for  (M<  =  100  GeV,  A///  =  100  GeV)  and  for 
{Mt  =  200  GeV,  Mjj  =  195  GeV)  axe  plotted.  These  two  couplings  are  the 
only  unknown  parameters  of  the  Standard  Model.  The  effects  of  changing  the 
values  of  Mt  and  Mff  in  the  analyses  of  the  running  of  the  other  parameters 
have  been  studied.  It  is  observed  that,  for  any  Mt  between  100  GeV  and  200 
GeV.  varying  Mj{,  while  maintaining  perturbativity  and  vacuum  stability,  did 
not  affect  appreciably  the  evolution  of  any  of  the  other  parameters.  However, 
changing  Mt  itself  showed  a  significant  difference  in  the  running  of  the  heav- 
ier quarks.  In  particular.  Fig.  10  shows  that  the  intersection  point  between 
the  bottom  quark  and  the  r  lepton  moves  down  to  a  lower  scale  for  higher 
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top  quark  masses.  This  is  expected  since  from  Eq.  (1.9)  one  can  see  that  the 
bottom  type  Yukawas  are  driven  down  by  an  increased  top  Yukawa.  This  is 
to  be  contrasted  with  the  SUSY  GUT  case  in  which  the  bottom  Yukawa  /3 
function  is  such  that  this  crossing  point  is  shifted  toward  a  higher  scale  with 
an  increased  top  mass.  The  relation  mj  =  m,-  (I)  is  the  most  natural  one  in 
the  SU{b)  theory  [44],  and  it  could  be  expected  to  be  valid  at  scales  where  the 
Standard  Model  gauge  couplings  axe  the  closest  to  one  another.  Its  validity 
is  examined  for  three  different  physical  values  of  the  top  and  Higgs  masses  in 
the  Standard  Model.  The  noteworthy  feature  of  Fig.  10  is  that  this  simplest 
of  the  SU{o)  relations  is  valid  at  an  energy  scale  many  orders  of  magnitude 
removed  from  that  at  which  the  gauge  couplings  tend  to  converge.  This  re- 
sult is  vastly  different  from  that  of  the  original  investigations  [45].  This  work 
improves  upon  that  work  by  including  two  loop  effects  in  the  running  of  the 
quark  Yukawas,  by  taking  into  account  the  full  Yukawa  sector,  and  most  im- 
portantly by  incorporating  QCD  corrections  in  the  extraction  of  the  bottom 
quark  mass.  Other  mass  relations  were  studied  in  the  SM  context  [46]:  (II) 
m^  =  3me,  3ms  =  m^  [47],  (III)  t^n6cabbibo  =  i^dl^s)^/'^  [48],  and  (IV) 
V^lj  =  {mc/mt)  '  [49].  There  is  no  scale  at  which  all  of  these  can  be  satisfied. 
The  scale  at  which  relation  (I)  tends  to  be  satisfied  does  not  coincide  with  that 
at  which  the  others  are  ^^alid.  Still  the  disagreement  is  never  too  large.  In  an 
SUi5)  SUSY  GUT  model,  the  equality  of  the  bottom  and  r  Yukawas  at  the 
scale  of  unification  will  be  used  to  get  bounds  on  the  top  and  Higgs  masses  [3]. 
Lastly,  the  running  of  the  CKM  angles  is  displayed  in  Fig.  11.  The  initial  data 
used  are  sin^i  =  0.2206,  sin  ^2  =  0.0298.  and  sin  ^3  =  0.0106.  Also  S  has  been 
taken  to  be  90°.  which  corresponds  to  the  case  of  maximal  CP  violation.   As 
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mentioned  in  Section  2.4,  the  evolution  curves  for  these  angles  are  effectively 
flat. 

In  the  present  case  of  the  Standard  Model,  it  is  found  that  two-loop  running 
of  the  parameters  does  at  times  improve  on  the  one-loop  running.  Indeed,  the 
differences  of  several  parameters  in  their  one-  versus  two-loop  values  at  various 
scales  have  been  tabulated,  for  the  cases  (Mt  =  100  GeV,  Mjj  =  100  GeV)  and 
{Mt  =  200  GeV,  Mfj  =  195  GeV).  Table  1  illustrates  the  difference  between 
one-loop  and  two-loop  running  in  the  ratio  mi^/mT,  for  the  three  scales  10 
GeV,  10  GeV,  and  10  GeV.  Clearly,  the  difference  between  one-  and  two- 
loop  results  is  more  pronounced  at  higher  scales,  as  expected.  Over  all  these 
scales  the  difference  is  never  less  than  5.5%.  Note  that  the  ratio  becomes  equal 
to  one  well  below  the  scale  of  grand  unification  as  noted  above  in  the  discussion 
of  Fig.  10.  Table  2  presents  a  similar  comparison  for  the  top  Yukawa.  Here, 
two  loops  represent  a  smaller  correction  with  the  difference  at  all  scales  always 
being  less  than  5%.  Finally,  Table  3  displays  the  same  analysis  for  a^  for 
the  case  M<  =  Mfj  =  100  GeV.  No  appreciable  deviation  from  the  tabulated 
values  is  observed  for  any  Mt  ^  200  GeV  (except  in  the  low  energy  regime 
where  the  difference  is  at  most  ~  4%).  At  scales  ^  M^,  the  inclusion  of  two 
loops  is  important  in  the  evolution  of  the  strong  coupling  (and  of  the  quark 
masses).  Indeed,  it  is  found  that  the  pure  QCD  three-loop  contribution  is  also 
significant  and  therefore  include  it  in  the  running  of  the  strong  coupling  and  of 
the  quark  masses  in  the  low  energy  region.  As  seen  in  this  table,  the  combined 
two  and  three  loops  in  the  low  energy  regime  account  for  a  17%  difference  at 
1  GeV  in  as. 

Although  in  the  cases  considered  in  these  last  two  tables  there  does  not 
appear  to  be  a  significant  difference  in  two-loop  over  one-loop  evolution  at 
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scales  above  Mz,  the  first  table  does  show  a  10%  difference  at  the  scale  10^° 
GeV.  The  effects  of  using  a  naive  step  approximation  vs.  a  proper  treatment  of 
thresholds  are  numerically  unimportant  for  the  cases  discussed  above.  Indeed 
they  are  less  important  than  the  two-loop  effects. 


Table  1:  mf^/mr 

Mt  =  100  GeV 

Mt  =  200  GeV 

102  GeV 

10'^  GeV 

10^6  GeV 

10^  GeV 

10^  GeV 

10^6  GeV 

One  loop 

1.879 

1.455 

0.8081 

1.868 

1.392 

0.6647 

Two  loop 

1.782 

1.348 

0.7336 

1.769 

1.285 

0.6047 

Table  2:  yt 

Mt  =  100  GeV 

Mt  =  200  GeV 

102  GeV 

10^  GeV 

10^6  GeV 

102  GeV 

lO'^  GeV 

10^^  GeV 

One  loop 

0.5405 

0.4160 

0.1928 

1.133 

0.9780 

0.7145 

Two  loop 

0.5405 

0.4071 

0.1842 

1.143 

0.9700 

0.6816 

Table  3:  Qs 

1  GeV 

102  GeV 

10^  GeV 

10^6  GeV 

One  loop 

0.3128 

0.1118 

0.07103 

0.02229 

Two 

and  three  loop 

0.3788 

0.1117 

0.07039 

0.02208 
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Figure  1.  Vacuum  stability  and  triviality  bounds  on  the  Higgs 
mass  for  Mt  =  100  GeV  giving  scales  of  expected  new 
physics  beyond  the  Standard  Model. 
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Figure  2.  Vacuum  stability  and  triviality  bounds  on  the  Higgs 
mass  for  Mt  =  125  GeV  giving  scales  of  expected  new 
physics  beyond  the  Standard  Model. 
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Figure  3.  Vacuum  stability  and  triviality  bounds  on  the  Higgs 
mELSS  for  Mt  =  150  GeV  giving  scales  of  expected  new 
physics  beyond  the  Standard  Model. 
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Figure  4.  Vacuum  stability  and  triviality  bounds  on  the  Higgs 
mass  for  Mt  =  200  GeV  giving  scales  of  expected  new 
physics  beyond  the  Standard  Model. 
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Figure  5.  Running  of  the  inverse  gauge  couplings  using  their  prop- 
agated experimental  errors  for  the  two-loop  case  only. 
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Figure  6.  Light  quark  and  lepton  masses  for  Mt  =  100  GeV  and 
Mh  =  100  GeV. 
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Figure  7.  Intermediate  quark  and  lepton  masses  for  Mt   =   100 
GeV  and  Mff  =  100  GeV. 
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Figure  8.  Heavy  quark  and  lepton  masses  for  Mt  =  100  GeV  and 
Mff  =  100  GeV. 
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Figure  9.  Top  Yukawa  and  scalar  quartic  couplings. 
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Figure  10.  Plot  of  mi,/mr  as  a  function  of  scale  in  the  Standard 
Model  for  various  top  and  Higgs  masses. 
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Figure  11.  CKM  mixing  angles  for  M<  =  100  GeV  and  M^  =  100 
GeV. 


CHAPTER  4 

THE  MSSM 


4.1  The  Supersvmmetric  Standard  Model 
In  the  minimal  supersvmmetric  extension  of  the  standard  model  (MSSM), 
every  particle  has  a  supersymmetric  partner,  their  spins  differing  by  a  half  [50]. 
Also  required  is  a  second  Higgs  field  with  opposite  hypercharge  to  the  first  as 
the  superpotential  cannot  contain  both  a  field  and  its  complex  conjugate.  The 
second  Higgs  is  also  needed  for  anomaly  cancellation  and  to  give  this  sector  a 
mass.  For  renormalizable  theories,  the  superpotential  can  have  at  most  degree 
three  interactions.  The  superpotential  for  the  MSSM  is  (suppressing  the  SU{2) 
and  Weyl  metrics) 

W  =  flYu^uQ  +  dYd^dQ  +  ^Ye^^L  +  jl^n^d  +  h.c.  ,  (4.1.1) 

where  the  hat  indicates  a  chiral  superfield  and  the  overline  denotes  a  left- 
handed  CP  conjugate  of  a  right-handed  field,  0  =  ia2tp*p-  The  usual  Yukawa 
interactions  are  accompanied  by  new  Yukawa  interactions  among  the  scalar 
quarks  and  leptons  and  the  Higgsinos  in  the  supersymmetric  Lagrangian.  There 
are  also  new  gauge  Yukawa  interactions  involving  the  gauginos.  The  new  purely 
scalar  interactions  form  the  scalar  potential  which  is  positive  definite  in  super- 
symmetric  theories.  The  scalar  potential  will  be  discussed  in  a  subsequent 
section. 

A  remarkable  aspect  of  supersymmetry  is  that  all  these  new  interactions 
require  no  new  couplings.  The  fi  term  is  only  present  to  avoid  a  Peccei-Quinn 
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(PQ)  symmetry.  Omitting  it  would  lead  to  exact  PQ  symmetry  and  to  a 
visible  axion  which  is  experimentally  ruled  out.  An  alternative  way  to  break 
the  PQ  symmetry  is  to  omit  the  jj,  term  and  add  an  explicit  soft  symmetry 
breaking  term,  m^f^u^d-  ^^^  A  can  be  interpreted  dynamically  as  essentially 
the  vacuum  expectation  value  of  a  singlet  chircd  superfield,  N,  through  the 
following  additional  interactions 

\N^r,^d  +  KNNN  .  (4.1.2) 

This  approach  also  provides  a  natural  explanation  for  jj.  ~  0{M\y).  The 
cubic  term  now  explicitly  breaks  the  PQ  symmetry.  Table  4  displays  the 
SU{3)  X  SU{2)  X  ^(1)  quantum  numbers  of  the  chiral  (ail  left-handed)  and 
vector  superfields  of  the  MSSM. 

4.2  Procedure 

Bounds  are  presented  for  the  mass  of  the  top  quark  in  a  minimal  supersym- 
metric  extension  of  the  Standard  Model  (MSSM)  with  minimal  Higgs  structure 
in  the  context  of  a  grand  unified  theory  (GUT)  by  numerically  evolving  the 
couplings  using  their  renormalization  group  equations.  This  analysis  improves 
on  previous  endeavors  by  taking  full  account  of  the  Yukawa  sector. 

In  the  expectation  that  the  Standard  Model  is  only  the  low  energy  mani- 
festation of  some  yet  unknown  GUT  or  of  a  possible  supersymmetric  (SUSY) 
extension  thereof,  the  three  couplings  ^3,  g2,  and  gi  corresponding  to  the  Stan- 
dard Model  gauge  groups,  SU(3)c  x  SU(2)i  x  ^(l)y-,  should  meet  at  some 
large  grand  unification  scale.  Using  the  accepted  values  and  associated  errors 
of  these  couplings  unification  is  observed  in  the  SUSY-GUT  case  but  not  in 
the  pure  GUT  case,  as  noted  by  several  groups  [1,51]  (see  Fig.  12).  However, 
this  should  not  be  viewed  as  proof  of  supersymmetry  since  given  the  \^lues  of 
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ai,  a2,  C13  at  some  scale,  and  three  unknowns  (the  value  of  a  at  the  unification 
scale,  the  unification  scale,  and  an  extra  scale  such  as  the  SUSY  scale)  there  is 
always  a  solution.  The  exciting  aspect  of  the  analysis  of  Ref.  1  is  the  numerical 
output,  namely  a  low  SUSY  scale,  Mgiigy  ■•  ^^^  ^  perturbative  solution  below 
the  Planck  scale  which  does  not  violate  proton  decay  bounds    [52]. 

Furthermore,  in  the  context  of  a  minimal  GUT  [44]  there  axe  constraints  on 
the  Yukawa  couplings  at  the  scale  of  unification.  One  first  restrict  oneself  to  an 
SU(5)  SUSY-GUT  [53]  where  yj  and  yr,  the  bottom  and  r  Yukawa  couplings, 
axe  equal  at  unification.  The  crossing  of  these  renormalization  group  flow  lines 
is  sensitive  to  the  physical  top  quark  mass,  A//.  This  can  be  seen  in  the  down- 
type  Yukawa  renormalization  group  equation  (above  Mgugy-i  for  example), 
from  which  the  evolution  of  yj  is  extracted,  since  the  top  contribution  is  large 
and  appears  already  at  one  loop  through  the  up-type  Yukawa  dependence 


dY 


-  =  J^Yd[  3Y/Yrf  +  Yu^Yu  +  TriZYj^Y^  +  YetVe} 


7  ifi  (4-2-1) 


dt         167r2 

15"'      -^z   ■     3 

where  Y^  ^^  g  are  the  matrices  of  Yukawa  couplings.  Demanding  that  their 
crossing  point  be  within  the  unification  region  determined  by  the  gauge  cou- 
plings allows  one  to  constrain  Mt.  This  yields  an  upper  and  lower  bound  for 
Mt  which  nevertheless  is  fairly  restrictive. 

There  is  a  threshold  at  MgjjgY  •  Here  the  matching  condition  is  the  naive 
one  of  simple  continuity  due  to  the  lack  of  knowledge  about  the  superparticle 
spectrum.  The  scale  is  taken  to  be  variable  to  account  for  this  ignorance. 

4.3  One  Light  Higgs  Limit 
The  simplest  implementation  of  supersymmetry  is  considered  and  the  cou- 
plings are  run  above  MsuSY   ^o  one  loop.     The  MSSM  is  assumed  above 
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^^SUSY-i  ^^d  ^  model  with  a  single  light  Higgs  scalar  below  it.  This  is  done  by- 
integrating  out  one  linear  combination  of  the  two  doublets  at  Mgugy  i  thereby 
leaving  the  orthogonal  combination  in  the  Standard  Model  regime  as  the  "Higgs 
doublet" 

^{SM)  =  *dCos^  +  ^«sin/3  ,  (4.3.1) 

where  $  =  ir2$*,  and  where  tan/3  is  also  the  ratio  of  the  two  vacuum  ex- 
pectation values  [vu/vj)  in  the  limit  under  consideration.  This  sets  boundary 
conditions  on  the  Yukawa  couplings  at  Mgugy .  Furthermore,  in  this  approx- 
imation the  quartic  self-coupling  of  the  surviving  Higgs  at  the  SUSY  scale  is 
given  by 

\{MsuSY)  =  \{9\  +  9h^os'^im  ■  (4.3.2) 

This  correlates  the  mixing  angle  with  the  quartic  coupling  and  thereby  gives 
a  value  for  the  physical  Higgs  mass,  Mfj.  Using  the  experimental  limits  on 
the  Mfj  further  constrains  some  of  the  results.  By  using  the  renormalization 
group  one  taJces  into  account  radiative  corrections  to  the  light  Higgs  mass  [54] 
and  hence  relax  the  tree  level  upper  bound,  Mf{  ~  M^    [55]. 

4.4  Initial  Data 
The  bounds  on  Mf  and  Mff  are  determined  by  probing  these  masses  de- 
pendence on  /3.  In  SUSY-5[^(5),  tan/?  is  constrained  to  be  larger  than  one  in 
the  one  light  Higgs  limit.  It  seems  natural  to  require  that  yt  >  y^  up  to  the 
uniJfication  scale  [56],  thereby  yielding  an  upper  bound  on  tan/3.  The  initial 
values  at  M^  for  the  gauge  couplings  are  taken  to  be    [1,57] 

ai  =  0.016887  ±0.000C40  , 

(12  =  0.03322  ±  0.00025  ,  (4.4.1) 

ninQ+0.004 
03  =  0.109_y^oQ5  , 
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where  GUT  normalization  for  ai  is  used.  The  following  set  of  four  quark 
running  masses  defined  at  1  GeV  by  the  Particle  Data  book  [21]  are  used 
mu  =  5.6  MeV,  m^  =  9.9  MeV,  rug  =  199  MeV,  and  rric  =  1.35  GeV.  For 
the  bottom  mass,  the  Gasser  and  Leutwyler  bottom  mass  value  of  5.3  GeV  at 
1  GeV  is  used  which  translates  into  a  physical  mass  of  A/j  =  4.6  GeV  [26]. 
To  probe  the  dependence  of  the  results  on  i\/j,  the  case  Mj  =  5  GeV  is  also 
studied,  this  is  the  typical  value  obtained  from  potential  model  fits  for  bottom 
quark  bound  states  [58].  The  effect  of  varying  Mgugy  is  also  investigated. 
Given  the  values  of  the  gauge  couplings,  unification  holds  for  SUSY  scales  up 
to  8.9  TeV  and  as  low  as  M\y.  For  empirical  reasons  solutions  below  My[/  were 
not  investigated. 

4.5  Analvsis  and  Results 
The  inclusion  of  supersymmetry  collapses  the  GUT  triangle.  This  is  il- 
lustrated in  Fig.  12  taking  Mgug}-  =  1  TeV.  As  mentioned  above,  a  range 
MiY  ^  Mgu^^  ^  9  TeV  will  achieve  unification  within  one  sigma  error.  From 
Fig.  13  (the  magnified  unification  region  of  Fig.  12),  one  determines  that  the 
lower  end  scale,  Mq^j,  of  the  unification  region  corresponds  to  an  aj,  value 
of  0.104  at  M^i  while  the  higher  end  scale,  Mqtjj,^  corresponds  to  a  value  of 
0.108  at  Mz  for  03.  It  is  found  that  the  unification  region  is  insensitive  to 
the  range  of  top,  bottom,  and  Higgs  masses  considered.  In  the  analysis  of  the 
bounds  for  Mt,  the  \'alues  for  qj  and  (12  ^^e  chosen  to  be  the  central  values 
since  their  associated  experimental  uncertainties  are  less  significant  than  for 
Q3.  Demanding  that  y^  and  j/r  cross  at  Mq^jj,  and  talcing  03  =  0.104  then 
sets  a  lower  bound  on  M^   Correspondingly,  demanding  that  y^  and  xjr  cross 
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at  MStjj.  and  talcing  03  =  0.108  yields  an  upper  bound  on  M<.  These  bounds 
are  found  for  each  possible  value  of  f3. 

Figure  14  shows  the  upper  and  lower  bound  curves  for  both  Mt  and  Mfj  as 
a  function  of  /?  and  for  Mgugy  —  1  TeV  and  Mj  =  4.6  GeV.  When  applicable 
the  current  experimental  limit  of  38  GeV  on  the  light  supersymmetric  neutral 
Higgs  mass  [59]  is  used  to  determine  the  lowest  possible  Mt  value  consistent 
with  the  model.  It  is  found  that  139  <  Mt  <  194  GeV  and  44  <  Mfj  <  120 
GeV.  The  sensitivity  of  these  results  on  Mgi/gy-  is  investigated  in  the  range, 
1.0  ±  0.5  TeV.  It  is  found  that  the  bounds  on  Mt  are  not  modified,  but  the 
upper  bound  on  the  Higgs  is  changed  to  125  GeV,  and  the  lower  bound  drops 
below  the  experimental  lower  bound. 

For  Mj,  =  5.0  GeV,  an  overall  decrease  in  the  top  and  Higgs  mass  bounds 
is  observed  116  <  Mt  <  181  GeV,  Mjj  <  111  GeV.  Varying  MsuSY  s-s  above 
modifies  the  respective  bounds.  The  top  mass  lower  and  upper  bounds  become 
113  and  119  GeV,  respectively.  The  upper  bound  on  Mjj  changes  to  115  GeV. 
The  results  of  the  analysis  are  displayed  for  the  extreme  case,  Mgjjgy  —  8.9 
TeV,  in  Fig.  15,  with  A/j  =  4.6  GeV.  This  only  significantly  changes  the  upper 
bound  on  M}j  to  144  GeV  compared  to  the  MgjjgY  —  1  TeV  case. 

yt  has  also  been  run  up  to  the  unification  region  and  compared  with  yj  and 
yr  to  see  what  the  angle  3  must  be  for  these  three  couplings  to  meet  [60],  as  in 
an  SO(IO)  or  E^  model  [61,62]  with  a  minimal  Higgs  structure.  It  is  clear  that 
this  angle  is  precisely  the  upper  bound  on  /?  as  described  earlier.  In  Fig.  16 
ytlyh  is  displayed  at  the  GUT  scale  as  a  function  of  tan/t?  for  Mgugy  =  1 
TeV  and  for  the  two  bottom  masses  considered.  If  one  demands  that  the  ratio 
be  one,  one  can  determine  the  mixing  angles  for  the  low  and  high  ends  of  the 
unification  region.  Then  going  back  to  Fig.  14.  one  finds  as  expected  a  much 
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tighter  bound  on  the  masses  of  the  top  and  of  the  Higgs.  Indeed,  for  M^  =  4.6 
GeV,  one  has  49.40  <  tan/^  <  54.98,  which  yields  162  <  A//  <  176  GeV  and 
106  <  Mh  <  111  GeV.  When  Mf,  =  5.0  GeV,  one  obtains  31.23  <  tanf3  < 
41.18,  which  gives  116  <  Mt  <  147  GeV  and  93  <  Mh  <  101  GeV. 

The  four  mass  relations  of  Section  3.7  were  also  studied  in  the  context  of 
the  MSSM  [46].  In  this  case,  these  relations  can  all  be  satisfied  at  the  scale  of 
gauge  coupling  unification.  However,  for  this  to  be  true,  several  things  must 
happen:  first  V^f,  must  be  larger  than  its  presently  measured  central  value  of 
0.043;  second  the  top  quark  mass  must  be  around  190  GeV  (if  it  is  lighter, 
then  agreement  dictates  that  l'f.5  should  be  larger  still);  third  the  Higgs  boson 
mass  should  be  around  120  GeV.  An  analysis  which  recently  appeared  in  the 
literature  has  reached  similar  conclusions    [63]. 

4.6  Comments 
To  improve  on  this  analysis,  one  should  implement  the  supersymmetric  two 
loop  beta  functions  and  the  corresponding  thresholds.  The  effects  of  soft  SUSY 
breaking  terms  should  be  investigated.  Also,  all  the  supersymmetric  particles 
have  been  integrated  out  at  the  same  scale.  It  would  be  interesting  to  study  the 
effect  of  lifting  this  restriction.  It  should  be  noted  that  the  bounds  on  the  top 
mass  are  very  similar  to  those  of  Ref.  56,  although  the  physics  is  very  different. 
These  issues  will  be  addressed  in  a  subsequent  chapter.  However,  given  the 
relative  crudeness  of  the  approximations  made  here,  it  is  remarkable  that  the 
experimental  bounds  on  the  p— parameter  were  satisfied.  A  sign  which  gives 
credence  to  the  program. 
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Table  4:  Particle  Quantum  Numbers  in  the  MSSM 
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Figure  12.  Plot  of  the  running  of  the  inverse  couplings.  The  dot- 
ted lines  above  and  below  the  solid  lines  represent  the 
experimental  error  for  each  coupling. 
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Figure  13.  The  plot  depicts  a  blow-up  of  the  area  around  the  unifi- 
cation point.  (Note  the  small  region  where  all  three  cou- 
plings intersect.  This  region  reduced  to  a  point  when 
MgugY  —  8-9  TeV  and  was  non-existent  above  that 
scale.) 
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Figure  14.  Plot  of  the  top  quark  mass,  M/,  and  of  the  Higgs  mass, 
Mfj,  as  a  function  of  the  mixing  angle  /9  for  the  highest 
value  of  a^  (high  curves)  and  the  lowest  value  of  0:3  (low 
curves)  consistent  with  unification  as  per  Fig.  12. 
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Figure  15.  Same  as  Fig.  14  for  MguSY  =  8-9  TeV  and  Mf,  =  4.6 
GeV. 
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Figure  16.  Plot  of  the  ratio  of  the  top  to  bottom  Yukawas,  yt/vb,  for 
two  different  bottom  masses  (soHd  and  dashed  curves) 
as  a  function  of  tan/3  for  the  highest  value  of  03  (high 
curves)  and  the  lowest  value  of  0:3  (low  curves)  consis- 
tent with  unification  as  per  Fig.  12. 


CHAPTERS 
SOFT  SYMMETRY  BREAKING 


5.1  Minimal  Low  Energy  Supergravitv  Model 
Since  no  super  particles  have  been  observed  experimentally,  supersymme- 
try,  if  truly  present  in  nature,  must  be  broken.  One  way  to  accomplish  this 
breaking  is  to  couple  the  standard  model  to  A'^  =  1  supergravity  (SUGRA). 
In  the  minimal  low  energy  supergravity  model  considered,  supersymmetry  is 
explicitly  broken  by  the  addition  of  supergravity  induced  soft  terms  (including 
gaugino  mass  terms) 

Vsoft  =  "^1^U«  +  "^L^i^d  +  Bfii^u^d  +  h.c.) 

i 

+  ^(  .'^dvt^lli^uQj  +  A'JY'Jdi^dQj  +  A^Yi^^i^iLj  +  h.c.  )  , 

«.; 

1    ^ 

Vgaugino  =  9  2-/  ^l^l^l  +  ^•<^-  ' 
"  /=1 

(5.1.1) 

where  Vgaugino  is  the  Majorana  mass  terms  for  the  gaugino  fields,  A/  (suppress- 
ing the  group  index),  corresponding  to  U{1),  5C/(2),  and  SU{3),  respectively. 
From  the  supersymmetry  algebra,  one  deduces  that  spontaneous  symmetry 
breaking  occurs  if  and  only  if  the  vacuum  energy  is  not  zero.  In  global  su- 
persymmetric  theories,  the  scalar  potential  is  a  sum  of  F  and  D  terms.  Su- 
persymmetry is  spontaneously  broken  if  either  the  vacuum  values  of  the  F 
term  (  [64])  or  D  term  (  [65])  are  non-zero.  A  consequence  of  the  spontaneous 
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symmetry  breaking  is  a  massless  fermion  in  analogy  with  the  breaking  of  an 
ordinary  global  symmetry.  We  will  assume  that  the  spontaneous  breaking  of 
the  local  N=l  supersymmetry  is  communicated  to  the  "visible"  sector  by  weak 
gravitational  interactions  from  some  "hidden"  sector.  Spontaneous  symmetry 
breaking  in  supergravity  occurs  via  the  super- Higgs  mechanism.  The  goldstone 
fermion,  or  goldstino,  associated  with  the  breaJcing  of  global  supersymmetry 
is  eaten  by  the  gravitino  thereby  providing  it  with  a  mass.  This  spontaneous 
symmetry  breaking  of  supergravity  manifests  itself  at  low  energy  as  explicit 
soft  breaking  terms  of  global  supersymmetry.  This  leads  to  a  common  (grav- 
itino) mass,  mg,  for  the  scalars  of  the  model  and  masses,  M/,  for  the  gauginos 
at  the  GUT  scale,  Mx-  By  assuming  gauge  coupling  unification,  we  can  take 
the  three  gaugino  masses  equal.  Furthermore,  the  trilinear  soft  couplings  Ay  , 
^y,  and  A''^  are  all  equal  to  a  common  value  .4o.  The  bilinear  soft  coupling  Bq 
may  be  related  to  A^  (Bq  =  A()  —  mg),  if  the  SUGRA  model  has  only  canonical 
kinetic  terms  for  the  chiral  superfields.  This  scenario  is  to  be  contrasted  with 
one  in  which  general  soft  breaking  terms  are  added  ad  hoc  to  the  Lagrangian. 
In  the  most  general  case,  there  axe  sixty-three  soft  symmetry  breaking  param- 
eters. Their  number  can  be  reduced  by  invoking  certain  symmetries  such  as 
flavor  and  family  blindness.  In  either  case,  all  of  these  couplings  will  evolve  to 
different  values  under  the  renormalization  group.  The  complete  scalar  poten- 
tial appears  as 

V^VF^-VD  +  V,,ft  ,  (5.1.2) 

where  Vp  contains  the  potential  contributions  from  the  F-terms 
Vf  =  \fiY,,Q  +  A$rf|2  +  i^^g  +  tY,Z  +  fi^ul^ 

+  \Y,Q^u\^  +  |Y^g$d  +  YeL^df  (5.1.3) 
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and  V/)  contains  the  potential  contributions  from  the  £)-terms 

2 

o 

(5.1.4) 
where  f  =  ( n  i '''2 1  ■''3 )  ^^^  the  St^(2)  Pauh  matrices  and  A  =  (Ai,...,A8)  are 
the  Gell-Mann  matrices.  In  general,  one  must  impose  constraints  on  the  pa- 
raxneters  to  avoid  charge  and  color  breaking  minima  in  the  scalar  potential. 
Some  necessary  constraints  have  been  formulated,  such  as 

Al<3{ml  +  ml  +  mlj  , 

Al<3{ml  +  ml  +  mlj,  (5.1.5) 

'4<3(mi  +  m2  +  m2j. 
However,  these  are  in  general  neither  sufficient  nor  indeed  always  necessary 
[66].  They  involve  very  specific  assumptions  about  the  spontaneous  symmetry 
breaking. 

5.2  Radiative  Electro- Weak  Breaking 
An  appealing  feature  of  the  models  being  considering  is  that  they  can 
lead  to  the  breaking  of  the  electro-weak  symmetry  radiatively.   The  one  loop 
potential  responsible  for  the  breaking  is 

Vl-loopif^)  =  Vtreeifi)  +  AFi(^)  ,  (5.2.1) 

where 

<4, 


647r2'  jj^ 


Al^l(/i)  -  — 2Str{A^^(ln-2-  -  -)} 
D47r^  jj.'^        2 

2  (5.2.2) 

=  ^E(-l)''^(2.p  +  l),nJ(ln5-|), 


P  ^ 
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where  A4    \s  the  field  dependent  squared  mass  matrix  of  the  model  and  irip  is 
the  eigenvalue  mass  of  the  p     particle  of  spin  Sp.  The  tree  level  potential  is 

^'2         t  t  -2        ^1         t  t  2  ^^-^-^^ 


where 


2  2,-2 


^2  =  "^!+/^^'  (5.2.4) 


Minimization  gives 


2        mj  —  m^  tan^  ,3 
2-'^"       tan2/3-l 


-I  =  "^LJr^:   ^  (5.2.5) 


where  m|  =  (y'"^  +  g2)v'^/2,  v^  =  v^  +  Vj^  and 


2 

where  tan/3  =  Vu/v^,  also 


5/i  =  -(mf  +  m^)  sin  2(3  ,  (5.2.6) 


-2  2   ,    ^AVi 


^i  =  "^i  +  -^  ,  (5.2.7) 


2 


where  t'l  =  u^,  U2  =  fu  and 

The  parameters  of  the  potential  are  taken  as  running  ones,  that  is,  they  vary 
with  scale  as  dictated  by  the  renormalization  group.  Because  the  one  loop 
correction  to  the  potential  is  not  negligible  at  all  scales,  an  appropriate  scale 
must  be  chosen  if  the  tree  level  formulas  are  to  be  valid  [67].  Indeed,  the  tree 
level  analysis  may  lead  to  incorrect  conclusions  about  the  regions  of  parameter 
space  that  yield  electroweak  breaking  and  consistent  scenarios.  When  AFi  is 
included,  the  value  of  ^  is  not  critical  as  long  as  /i  is  in  the  neighborhood  of 
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Although,  as  stated,  the  tree  level  results  cannot  always  be  trusted,  one 
c£Ln  get  some  idea  under  what  conditions  electro-weak  breaking  occurs.  The 
renormalization  group  evolution  of  mi  (see  Appendix  C)  can  be  such  that  it 
turns  negative  at  low  energies  if  the  top  quark  mass  is  large  enough,  whereas 
m?    runs  positive.   From  Eq.  (5.2.3),  the  scale  at  which  this  occurs  is  set  by 


the  condition 


mlM  mlM  -  ml(fih)  =  0  .  (5.2.9) 


If  the  free  parameters  are  adjusted  properly,  then  the  correct  value  of  the 
Z^  mass  {Afz  —  91.17  GeV)  can  be  achieved.  At  tree  level  there  is  another 
critical  scale  that  must  be  considered.  Again  from  Eq.  (5.2.3),  it  is  evident  that 
the  potential  becomes  unbounded  from  below  along  the  equal  field  (neutral 
components)  direction,  if 

m\{^is)  +  ml{yis)  <  2m^(^s)  .  (5.2.10) 

Since  rri^m^  —  m^  >  0  impHes  m^  -|-  m2  >  27713,  condition  ((5.2.10))  can  only 
occur  at  scales  lower  than  condition  ((5.2.9)),  so  ^s  <  /;^.  As  mentioned  above, 
when  working  with  the  tree  level  potential  an  appropriate  scale  to  minimize  it 
must  be  used  at  which  one  loop  corrections  may  be  safely  neglected.  Gamberini 
et  al.     [67]  give  a  prescription  for  the  choice  of  this  scale. 

In  the  present  work,  the  one  loop  corrections  are  incorporated.  Contribu- 
tions from  the  third  generation  are  included,  that  is,  the  top-stop,  bottom- 
sbottom,  and  tau-stau.  The  one  loop  effective  potential  should  be  constant 
against  the  renormalization  group  to  this  order.  The  Z^  mass  is  chosen  as  the 
scale  at  which  to  evaluate  the  minimization  conditions  Eqs.  (5.2.9)-(5.2.10). 
Equation  (5.2.5)  can  be  written 

o  ^A    —  ^\   tan^  8       1 
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where  mj^  ^  =  ^l^  ^  +  d^Vi/dvl^^  and  used  to  solve  for  il{Mz)  (see  Numerical 
Procedure  Section)  given  the  value  of  all  the  relevant  parameters  at  Mz-  A 
choice  for  the  sign  of  /i  must  be  made  [jj.  is  multiplicatively  renormalized;  see 
Appendix  C).  It  is  evident  that  no  adequate  minimum  exists  if  the  parameters 
are  such  that  fi    <  0. 

5.3  Sfermion  Masses 
The  mass  matrices  for  scalax  matter  are  constructed  from  Eq.  (5.1.2).  For 
example,  in  the  up  squark  sector  the  relevant  mass  matrix  appears  as 

\      RiLj  LiRj  / 

where  i,j  =  1,2,3  are  flavor  indices  and 

^II,  =  ^Q.^ij  +  vl(YiY,)ij  -  ^{vj  -  vl)(YiuL)g''^  -  T^{uL)gl)6i^  , 
^luR,  =  rr^i^ij  +  ^'(ySyJ,-,-  -  \{vl  -  vl){Y{uR)g''^)8i^  , 
Mr,Lj  =  i^^dYi^  +  VuA^Yi^  , 

m"^       =  A/2* 
Iji  Hj  Rj  Li 

(5.3.2) 
Note  in  Table  4  that  Y  =  Q  -  T^  in  the  adopted  notation.  Similar  matrices 
follow  for  the  other  sfermions.  These  mass  formulas  as  well  as  the  ones  to 
follow  are  given  in  terms  of  running  parameters.  The  domain  of  validity  of 
these  formulas  is  at  low  energies  (~  Mz)  with  the  parameters  taking  on  their 
renormalization  group  evolved  values  at  this  scale. 

5.4  Higgs  Masses 
If  the  following  notation  is  employed 

-M)  ■-(:!)■ 
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then  the  physical  masses  of  the  Higgs  at  tree  level  are  calculated  from  the 
following  three  matrices 

2a(3J<^?)5W5)       2  ^    "         ^     ^^  (5  4  1) 

,1     2     •    o«  /cot/?        -1    \ 
+  -m^sm2^(^    -1       tan/? J 

^^^<ree  1  ,  ,2        ■     ^nli^xv^  1      \ 
x~  =  x-Mu    sm2/3         1  .  o 

where  A/|  =  mj  +  m^,  A/^^  =  MJ  +  m^,  and  mjy  =  g\v^  12.  The  eigenval- 
ues for  the  first  matrix  are  0  corresponding  to  the  Goldstone  boson  and  M\ 
corresponding  to  the  CP  odd  scalar.  The  second  matrix  gives  the  masses  of 
the  light  and  heavy  Higgs  bosons 

^EM  =  \k'^\  +  rn\)  ±  ^(A/2  +  m|)2  -  m\m\  cos2  2/?]  .  (5.4.2) 

This  tree  level  result  predicts  A//,  <  Mz-    One  loop  calculations  show  that 
this  need  not  be  the  case    [68].  The  third  matrix  has  eigenvalues  0  and  M\t 
corresponding  to  a  massless,  charged  Goldstone  boson  and  a  charged  scalar.  In- 
cluding Eq.  (5.2.2)  in  the  calculations  leads  to  corresponding  one  loop  versions 
of  these  masses    [69,70]. 

5.5  Chargino  Masses 
The  following  four  terms  contribute  to  the  chargino  masses 


-i^g2^V-l^~^nW'  -  i^flg2^'^^-~^  iW'  -  li^u^d  +  \miW'W'  +  h.c.  .   (5.5.1) 
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The  first  two  terms  are  the  supersymmetric  Yukawa-gauge  terms.     Letting 
A^  =  {W2  ±  iWi)/\/2,  the  mass  matrix  follows 


(a+     ^+     A-     4>^)\ 


/     0           0  M2  -g2Vd\ 

0            0  g2Vu  -M 

M2  g2Vu        0  0 

\-g2^'d      -A         0  0      / 


^-^\ 


Diagonalization  yields  two  charged  Dirac  fermions  with  masses 


±  y^(M|  +  /i2  +  2m2^)2  -  4(M2/i  -  m'J^  sin  2/3)2]  . 


(5.5.2) 


(5.5.3) 


5.6  Neutralino  Masses 
Contributing  to  the  neutralino  masses  are  the  terms  in  Eq.  (5.5.1)  and 


-z9'V2^U+\)^uB  -  Vv^$^(-^)M  +  ^MiBB  . 


(5.6.1) 


The  neutralino  mass  matrix  follows 

f  -A/i 
1 


(zB     iW^     <^0     ,^0)- 


0 


-M2 


g  Vd       _21M 


V2  V2 

g2t'd  ■?2t'ti 

0 


/^ 


0 


/ 


(5.6.2) 


5.7  Boundary  Conditions  at  A/y 
In  this  work,  the  modified  minimal  subtraction  scheme  (MS)  of  renormal- 
ization  is  employed.  The  parameters  of  the  Lagrangian  are  not  in  general  equal 
to  any  corresponding  physical  constant.  For  example,  in  the  case  of  masses, 
except  for  those  of  the  bottom  and  top  quark  (see  Eq.  (3.4.1)),  all  other  physi- 
cal masses  will  be  determined  from  their  corresponding  running  masses  by  the 
simpler  relation 


M  =  rn(^) 


H=M 


(5.7.1) 


57 

The  renormaiization  group  f3  functions  of  the  gauge  and  Yukawa  coupUngs 
have  been  calculated  to  two  loops  without  making  any  approximations  in  the 
Yukawa  sector  of  the  model.  These  have  been  included  in  Appendix  C. 

Because  SUGRA  models  make  certain  simplifying  predictions  about  the 
soft  parameters  at  the  unification  (Planck)  scale,  the  evolution  of  the  renor- 
maiization group  equations  is  initiated  at  this  scale.  It  has  been  demonstrated 
that  the  introduction  of  supersymmetry  leads  to  gauge  coupling  unification  at 
approximately  ~  10      GeV.  Therefore  one  takes  Mx  =  10      GeV. 

At  the  unification  scale,  Mx,  all  scalars  have  a  common  mass 

^Q,{Mx)  =  rriuiiMx)  =  rn^^iMx)  =  mi^{Mx) 

(5.7.2) 
=  metiMx)  =  m^^{Mx)  =  m^,{Mx)  =  rriQ  , 

as  do  the  gauginos 

MiiMx)  =  M2iMx)  =  MziMx)  =  mi  .  (5.7.3) 

2 

The  trilineax  soft  scalar  couplings  will  be  talcen  diagonal  and  equal  at  Mx 

A^iiMx)  =  A'jiMx)  =  .^J{Mx)  =  -4o  .  (5.7.4) 

Also  the  bilinear  soft  scalar  coupling  and  the  mixing  mass  at  A/y  are  given  by 

B{Mx )  =  5o,  fi{Mx)  =  fiQ  .  (5.7.5) 

Furthermore,  to  constrain  the  parameter  space,  the  bottom  and  tau  masses 
will  be  assumed  equal  at  Mx 

mbiMx)  =  rnr{Mx)  ■  (5.7.6) 

Also  for  the  sake  of  CPU  time,  we  will  take  the  masses  of  the  two  lightest 
families  to  be  zero  in  all  the  runs. 
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5.8  Numerical  Procedure 
There  are  seven  free  parameters  in  the  model  considered.  These  are  .4o, 
Bq,  rriQ,  "^1/25  AOi  tan /^,  and  mt.  The  two  minimization  constraints  (5.2.5) 
and  (5.2.6)  reduce  this  set  to  five,  which  are  taken  to  be  Aq,  mg,  rn\i2i  tan/3, 
and  m^  In  the  present  framework,  jBq  and  fiQ  will  be  determined  using  the 
numerical  solutions  routines  described  in  Appendix  E  in  conjunction  with  the 
minimization  of  the  one  loop  effective  potential  at  Mz-  Minimization  at  Mz 
will  give  B{Mz)  and  fi[Mz)-  To  arrive  at  Bq  and  fiQ  (their  corresponding 
values  at  My),  the  solution  routine  are  employed  as  follows.  A  guess  for  Bq 
and  jj,Q  is  made  at  A/y  and  then  the  parameters  of  the  model  are  run  to  Mz 
at  which  scale  the  evolved  value  of  B  is  compared  to  the  minimization  out- 
put value  for  B  at  M^.  The  saxne  is  done  for  fi.  If  the  compared  values 
agree  to  some  set  accuracy,  then  Bq  and  fiQ  are  the  required  values.  Other 
analyses  that  also  extract  B{Mz)  and  fi{Mz)  simply  evolve  these  two  param- 
eters via  their  renormalization  group  equations  back  to  Mx  to  find  Bq  and 
fiQ  relying  on  their  near  decoupling  from  the  full  set  of  renormalization  group 
equations.  Notice  from  Eq.  (5.2.11)  that  a  choice  must  be  made  for  the  sign 
of  ft.  To  constrain  the  parameter  space  further,  the  bottom  quark  and  tau 
lepton  masses  will  be  taken  equal  at  Mx-  This  equality  is  a  characteristic  of 
many  SUSY-GUTs.  This  constrains  the  model  to  four  free  parameters,  .-lo, 
mo,  rni/2i  and  tan/3.  Demanding  that  mj(MY)  =  rnT{Mx)  and  achieving  the 
correct  physical  masses  for  the  bottom  quark  and  tau  lepton  fixes  the  mass  of 
the  top  quark  which  affects  the  evolution  of  the  bottom  Yukawa  significantly. 
Gauge  coupling  unification  shall  be  assumed,  an  assumption  which  appears 
reasonable  when  one  considers  SUSY  models  with  SUSY  breaking  scales  ^  10 
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TeV.  The  solution  routines  could  be  used  to  find  the  precise  (and  similax)  val- 
ues of  ai,  a2i  ^^^  <^3  ^^  -^^.Y  ^^^^  '^^^^^  evolve  to  the  experimentally  known 
values  at  M^,  however  this  increases  the  CPU  time  considerably.  Therefore 
some  precision  shall  be  sacrificed  in  their  Mz  values  by  taking  them  exactly 
equal  at  My-  This  is  already  a  theoretical  oversimplification  since  one  does 
not  expect  the  gauge  couplings  to  be  exactly  equal  due  to  threshold  effects 
at  the  GUT  scale.  It  is  found  that  for  all  ccises  studied,  the  common  value 
a'^^iMx)  =  <^2^^(-^-^.y)  =  ^s'H-^x)  =  25.31  leads  to  errors  no  bigger  than 
1%,  5%,  and  10%  in  ai(M^),  a2{Mz),  and  a:i{Mz),  respectively.  This  is  not 
so  bad  considering  that  the  (combined  experimental  and  theoretical)  errors  on 
a^{Mz)  from  some  processes  can  be  as  large  as  this. 

It  is  well  known  that  there  is  a  fine  tuning  problem  inherent  in  the  ra- 
diatively  induced  electro-weak  models.  For  certain  values  of  the  parameters, 
the  top  quark  mass  must  be  tuned  to  an  "unnaturally"  high  degree  of  accu- 
racy to  achieve  the  correct  value  of  Mz-  This  problem  is  generally  handled 
by  rejecting  models  that  require  "too  much"  tuning.  The  amount  of  tuning  is 
usually  defined  quite  arbitrarily.  The  usual  procedure  is  to  define  fine  tuning 
parameters 

--'— ^ ^'  (5.8.1) 


'M|  dx] 


where  ii  are  parameters  of  the  theory  such  as  mg,  rn\/2i  f^>  o^  '^^  One  then 
demands  that  the  c,  be  less  than  some  chosen  value  that  is  typically  taken  to  be 
10.  The  differences  in  using  the  tree  level  vs.  one  loop  effective  potential  were 
analyzed  to  some  extent.  The  results  agree  generally  with  those  of  Ref.  67. 
Moreover  it  is  found  that  the  fine  tuning  problem  is  exacerbated  in  the  tree 
level  analysis.  The  basis  for  the  '"theoretical"  fine  tuning  problem  can  be  seen, 
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if  one  makes  some  simplifying  assumptions,  in  the  dependence  of  M^r  on  the 
top  quark  Yukawa  coupUng  yt    [71] 

Mw  ~  Mxe'^/y'^  .  (5.8.2) 

In  the  tree  level  analysis,  one  encounters  another  fine  tuning  problem.  The 
vacuum  expectation  value  coming  from  the  minimization  conditions  changes 
rapidly  from  0  to  infinity  over  the  interval  (;is,/X{,).  Using  the  prescription  of 
Ref.  67  for  the  scale  (i  at  which  to  adequately  minimize  the  tree  level  potential, 
find  v{fi),  and  thereby  arrive  at  a  value  for  M^,  one  finds  that  although  a  small 
variation  in  yt{Mx)  may  lead  to  a  small  variation  in  ^j,  the  steepness  in  the 
tree  level  vacuum  expectation  value  will  lead  to  a  large  variation  in  the  value 
of  v{jl)  and  therefore  in  Mz-  Hence,  in  the  tree  level  analysis,  solutions  which 
may  be  within  the  bounds  of  the  "theoretical"  fine  tuning  may  nevertheless 
display  a  fine  tuning  aspect  because  of  this  "tree  level"  fine  tuning  problem. 
The  use  of  the  one  loop  effective  potential  levels  off  the  vev  around  the  M^ 
scale.  The  vevs  depend  on  scale  through  wave  function  renormalization  effects 
which  are  never  large  as  can  be  seen  from  the  form  of  the  renormalization 
group  equations  for  the  vevs  in  Appendix  C. 

In  the  present  framework,  solutions  shall  also  be  rejected  based  on  fine 
tuning  considerations,  however  the  present  method  differs  somewhat  from  the 
usual  one  in  that  it  is  incorporated  in  the  solution  routine  described  above. 
The  routine  is  an  iterative  one  which  determines  the  convergence  properties 
of  the  solution  which  reflect  an  inherent  fine  tuning.  If  the  convergence  is  too 
slow,  the  solution  will  be  rejected.  Effectively  any  solution  which  the  computer 
cannot  pinpoint  within  an  allotted  number  of  iterations  is  rejected.  Given 
values  for  .-lo,  mg,  '^i/2)  tan/i,  and  sign(/i),  the  solution  routines  search  for 
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the  values  of  v{Mx),  mi,T.{Mx),  mt(Mx},  Bq,  and  JIq.  The  process  by  which 
Bq  and  fiQ  are  found  was  described  above.  The  remaining  three  parameters 
are  determined  similarly.  The  routine  makes  a  guess  for  v(My),  m5^(MY), 
and  mt(Mx),  then  the  full  renormalization  group  equations  are  evolved  to  1 
GeV  calculating  superparticle  threshold  masses  in  the  process  and  minimizing 
the  one  loop  effective  potential  at  Mz-  The  merits  of  the  guess  for  v{Mx), 
mi  t{Mx),  and  mt(Mx)  is  assessed  by  comparing  the  resulting  values  of  A/^, 
mr(l  GeV),  and  mft(l  GeV)  with  the  expected  ones.  The  process  is  iterated 
until  the  correct  values  axe  achieved  to  within  some  tolerance. 

5.9  Thresholds 
In  the  minimal  low  energy  supergravity  model  being  considered,  the  super 
particle  spectrum  is  no  longer  degenerate  as  in  the  simple  global  supersymme- 
try  model  in  which  all  the  super  particles  are  given  a  common  mass.  Mgugy  • 
In  the  simple  case,  one  makes  one  course  correction  in  the  renormalization 
group  evolution  at  MgjjgY-  I"  ^-he  model  with  soft  symmetry  breaking,  the 
nondegenerate  spectrum  should  lead  to  various  course  corrections  at  the  super 
particle  mass  thresholds.  To  this  end,  the  renormalization  group  /3  functions 
must  be  cast  in  a  new  form  which  makes  the  implementation  of  the  thresholds 
effects  (albeit  naive)  evident  (see  Appendix  B).  Since  the  MS  renormalization 
group  equations  are  mass  independent,  particle  thresholds  must  be  handled  us- 
ing the  decoupling  theorem  [39],  and  each  super  particle  mass  has  associated 
with  it  a  boundary  between  two  effective  theories.  Above  a  particular  mass 
threshold  the  associated  particle  is  present  in  the  effective  theory,  below  the 
threshold  the  particle  is  absent. 
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The  simplest  way  to  incorporate  this  is  to  (naively)  treat  the  thresholds  as 
steps  in  the  particle  content  of  the  renormalization  group  (3  functions.  This 
method  is  not  always  entirely  adequate.  For  example,  in  the  case  of  the  SU{2) 
gauge  coupling  there  will  be  scales  in  the  integration  process  at  which  there  are 
effectively  a  half  integer  number  of  doublets  using  this  method.  Nevertheless, 
this  method  should  yield  the  correct,  general  behavior  of  the  evolution.  It  is 
a  simple  means  of  implementing  the  smeziring  effects  of  the  non- degenerate 
super  particle  spectrum.  The  determination  of  the  spectrum  of  masses  is  done 
without  iteration  as  is  common  in  other  analyses.  The  method  employed  in 
this  work  deduces  the  physical  masses  by  solving  the  equation  m(^)  =  jj,  for 
each  superparticle  in  the  process  of  evolving  from  Af  y  •'^  1  GeV.  The  usual 
iterative  method  requires  several  runs  to  find  a  consistent  solution. 

5.9  Analvsis  and  Results 

The  tremendous  computing  task  involved  in  analyzing  the  full  parameter 
space  of  the  soft  symmetry  breaking  models,  using  the  methods  described 
as  designed,  would  be  far  too  time  consuming  given  the  available  computing 
facilities.  Therefore,  in  the  following  analysis,  some  simplifications  will  be  made 
in  the  procedural  method.  First,  only  the  heaviest  family  of  quarks  and  leptons 
will  have  non- zero  mass.  Second,  as  stated  previously,  the  value  of  the  strong 
coupling  at  M^  will  be  allowed  to  vary  from  its  central  value  of  .113  by  at  most 
10%.  This  translates  into  a  similar  error  in  the  bottom  quark  mass.  Third, 
the  allotted  number  of  Runge-Kutta  steps,  involved  in  numerically  integrating 
the  renormalization  group  equations,  will  be  cut  down  to  ~  100. 

The  present  analysis  will  restrict  itself  to  a  subclass  of  soft  symmetry  break- 
ing models  that  have  two  or  three  soft  parameters  equal  to  zero  at  A/y-  One 
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such  class  follows  from  the  no-scale  model  [72]  and  has  Aq  =  tuq  =  0.  The 
strict  no-scale  model  has  .4o  =  mo  =  Bq  =  0.  Another  class  coming  from 
string  derived  models  has  Aq  —  Bq  =  0.  Since  Bq  is  not  a  free  parameter  in 
the  procedure  adopted  in  this  work  (see  Section  5.8),  Bq  =  0  results  must  be 
inferred. 

Because  the  GUT  inspired  constraint,  mi,{Mx)  =  mriMx),  is  enforced 
in  this  analysis,  the  results  will  depend  on  the  mass  of  the  bottom  quark. 
Most  results  will  be  reported  for  the  case  ^^(1  GeV)  =  6.00  GeV,  but  lower 
mass  (5.70  GeV)  and  higher  mass  (6.33  GeV)  cases  were  also  studied.  The 
running  value  of  6.00  GeV  for  mj(l  GeV)  corresponds  to  a  physical  bottom 
mass  Mj  =  4.85  ±  .15  GeV,  with  the  uncertainty  coming  from  the  error  in  the 
strong  coupling,  as  discussed  above.  The  results  of  Chapter  4,  indicate  that 
the  top  quark  mass  predictions  increase  with  decreasing  bottom  mass,  and  the 
present  analysis  corroborates  this. 

Given  the  no-scale  or  string-inspired  cases,  the  phase  space  is  explored  by 
setting  .4o  =  ttiq  =  0  ov  Aq  =  0,  respectively,  and  coarse  graining  the  remaining 
hyperslice.  tan/?  was  most  commonly  coarse  grained  as  2,  5,  and  10.  The  mass 
of  the  lowest  supersymmetric  particle  (LSP),  when  it  is  a  neutralino,  is  observed 
to  be  correlated  with  the  value  of  mw2.  Cosmological  considerations  indicate 
that  the  LSP  must  be  neutral  and  colorless  and  have  a  mass  less  than  ~  200 
GeV.  Points  in  parameter  space  that  lead  to  LSPs  other  than  neutralinos  with 
masses  less  than  ~  200  GeV  (or  sneutrinos)  are  rejected. 

Figure  17  shows  a  slice  of  the  available  pha^e  space,  in  the  no-scale  case, 
plotted  against  the  top  quark  mass.  The  present  unofficial  bound  on  the  top 
quark  mass  is  ~  120  GeV.  The  figure  indicates  that  the  top  quark  cannot  have 
a  mass  greater  than  ~  132  GeV  in  this  model  if  7775(1  GeV)  =  6.00  GeV 
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(This  upper  bound  is  raised  to  ~  160  GeV,  if  mi,(l  GeV)  =  5.70  GeV,  and  the 
model  is  ruled  out,  if  mf,(l  GeV)  =  6.33  GeV).  The  figure  also  indicates  that 
190  GeV  ^  rni/2  ^  265  GeV;  however,  nothing  significant  can  be  said  about 
the  value  of  tan/3.  The  points  to  the  right  of  the  allowed  region  were  found 
to  lead  to  charged  LSPs  (fn)  and  are  not  displayed.  Points  to  the  left  of  the 
region  do  not  lead  to  electro-weak  breaking.  Furthermore,  as  Figs.  14  and  15 
indicated,  the  top  quark  mass  does  not  grow  monotonically  with  tan/?,  rather 
it  reaches  a  maximum  for  some  high  tan  /3  value  then  decreases.  For  the  choice 
sign(^)  =  -,  the  top  quark  mass  upper  bound  falls  below  the  experimental 
limit,  so  this  case  is  ruled  out. 

Analysis  of  the  available  data  indicates  that  the  strict  no-scale  case  lowers 
the  top  quark  mass  upper  bound  shghtly  to  ~  128.  However,  there  is  an  allowed 
range  of  3  -  9  for  tan/3.  Figure  IS  is  a  plot  of  the  evolution  of  the  left-handed 
and  right-handed  stop,  sbottom,  stau,  and  gaugino  soft  mass  parameters  for 
the  particular  strict  no-scale  case  mi/2  =  240  GeV,  tan/3  =  8.3,  sign(//)  -  +. 
The  resulting  spectrum  of  super  particle  masses  is  presented  in  Table  5.  The 
LSP  is  a  neutralino  in  this  case  and  has  a  mass  of  92  GeV.  The  top  quark  mass 
is  just  above  the  experimental  limit  at  126  GeV,  and  the  Higgs  boson  mass  is 
77  GeV. 

The  strict  no-scale  case  is  a  special  case  of  the  string  inspired  one.  therefore 
the  128  GeV  top  quark  mass  upper  bound  is  not  expected  to  decrease  but  rather 
to  increase  in  the  "stringy"  case.  The  data  indicate  that  top  quark  masses  as 
high  as  ~  150  GeV  are  possible.  Finally,  a  representative  stringy  scenario  is 
presented.  Figure  19  displays  the  rariation  of  Mf  and  Bq  with  tan  3  for  the 
slice  of  parameter  space  with  .4o  =  0,  7?7o  =  100  GeV,  and  mj/2  =  250  GeV. 
The  scenario  inferred  from  the  figure  has  Mf  =  130  and  tan/3  =  7.8. 
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Table  5 

-4o 

0 

mo 

0 

"^1/2 

240 

Ao 

171 

tan/3 

8.3 

sign(^) 

+ 

Mt 

126 

d 

473  -  506 

u 
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Figure  17.  Slice  of  parameter  space  in  no-scale  case  displaying  A/< 
against  the  common  gaugino  mass  mw2-  The  points 
represent  radiative  electro-weak  breaking  solutions  with 
rrif,  =  tut  at  Mx  and  a  neutralino  LSP  with  mass  ^  200 
GeV. 
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Figure  18.  Evolution  of  squark  (dashes),  slepton  (dots),  and  gaug- 
ino  (solid)  soft  mass  terms  for  .4o  =  0,  mg  =  0,  mw2  = 
240  GeV,  and  tan /^  =  8.3. 
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Figure  19.  Mt  (solid)  and  Bq  (dashes)  plotted  against  tan /?  for 
.4o  =  0,  mo  =  100  GeV,  and  m^/j  =  250  GeV  to  infer 
the  Bq  =  0  case. 


CHAPTER  6 
CONCLUSIONS 

This  work  has  studied,  using  one  and  two  loop  renormahzation  group  l3 
functions,  both  the  standard  model  and  its  minimal  supersymmetric  extension. 
The  parameters  of  these  models,  i.e.,  the  gauge  couplings,  the  quark,  squark, 
lepton,  and  slepton  masses,  the  Yukawa  sector  mixing  angles  and  phase,  the 
scalar  quartic  coupling,  and  soft  symmetry  breaking  parameters,  were  run  over 
scales  ranging  from  1  GeV  to  Planck  mass.  The  aspects  of  the  standard  model 
and  its  minimal  supersymmetric  extension  were  reviewed. 

In  the  standard  model  case,  thresholds  effects  were  served  well  by  naive  step 
functions.  The  more  sophisticated  implementation  using  one  loop  matching 
functions  did  not  significantly  improve  on  the  step  function  method.  The 
difference  between  two  loops  and  one  loop  represented  a  more  significant  effect, 
albeit  sometimes  small.  Plots  exhibiting  these  different  features  for  all  the 
parameters  of  the  Standard  Model  were  included.  Gauge  coupling  unification 
was  shown  to  fail,  and  some  interesting  mass  and  mixing  angle  relations  were 
considered,  but  failed  to  hold  simultaneously  at  a  common  scale. 

In  the  supersymmetric  case,  in  which  gauge  coupling  unification  is  achieved, 
bounds  on  the  top  quark  and  Higgs  boson  masses  were  determined  using  the 
SU{5)  inspired  constraint  that  the  bottom  and  tau  meisses  be  equal  at  the 
scale  of  unification.  Remarkably,  the  top  quark  bounds  were  consistent  with 
the  p  parameter  bounds.  In  supersymmetric  context,  the  mass  relations  could 
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cdl  be  accommodated  at  the  scale  of  unification  provided  the  top  quark  mass 
was  high  ~  190. 

Minimal  low  energy  supergravity  models  were  considered.  They  have  the 
appealing  feature  that  the  electro-weal;  symmetry  is  radiatively  broken  for  cer- 
tain ranges  of  the  soft  breaking  parameters.  The  study  of  specific  models,  with 
some  soft  parameter  fixed,  resulted  in  upper  bounds  for  the  top  quark  mass. 
No-scale  models  in  which  only  gaugino  masses  provide  global  supersymmetry 
breaking  yield  top  qucurks  with  masses  less  thaxi  ~  130  GeV.  The  results  are 
sensitive  to  the  value  of  the  bottom  quark  mass.  Lower  bottom  quark  masses, 
within  the  experimental  uncertainty,  lead  to  higher  top  quark  upper  bounds. 
In  these  models,  the  ratio  of  vacuum  expectation  values  of  the  two  Higgs  fields 
is  expected  to  be  larger  than  ~  72°. 


APPENDIX  A 
THE  STANDARD  MODEL  i3  FUNCTIONS 

In  this  appendix,  the  renormalization  group  /3  functions  of  the  Standard 
Model  are  compiled  [6].  These  have  appeared  in  one  form  or  another  in  various 
sources.  Effort  has  gone  into  confirming  their  validity  through  a  comparative 
analysis  of  the  literature.  The  main  source  is  Ref.  7.  Following  their  conven- 
tions. 

£  =  uYu^^Q  +  dV^^^Q  +  eYe$^I  +  h.c.  -  ^A($t$)2  ,  (A.l) 

where  flavor  indices  have  been  suppressed,  and  where  Q  and  L  are  the  left- 
handed  quark  and  lepton  SU(2)  doublets,  respectively. 


$  and  $  are  the  Higgs  scalar  doublet  and  its  SU(2)  conjugate 

•^  =  (  lo  )    ,  $  =  zr2$  .  (A.3) 


Ur,  da,  and  Cr  are  the  quark  and  lepton  SU(2)  singlets,  and  Y^^^Le  are  the 
matrices  of  the  up-type,  down-type,  and  lepton-type  Yukawa  coupHngs. 
The  /3  functions  for  the  gauge  couplings  are 

dt  "'l67r2       Z.''^-'(i67r2)2 

(A.4) 
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where  t  —  In/i  and  /  =  1,  2,  3,  corresponding  to  the  gauge  group  SU{Z)c  x 
SU{2)i  X  U{\)y  of  the  Standard  Model.  The  various  coefficients  are  defined 
to  be  ■  ^ 


Vo 


and 


136 

T 


6i  =  --rin  - 


10  ' 


22      4 

1 

i'2  =  y  -  3", 

~  6 

4 
63  =  11  -  -ng 

0   \ 
102/ 

-   Ug 

3       49 
5      T 
44        . 
\T5      4 

11  \ 

(  9 

17      1       3 
,   T(J     5     5   , 
'Q/)  =14      4      A   1    '  with/  =  rt  ,  c?  ,  e  , 


(A.5) 


(A.6) 


(A.7) 


with  n. 


5"//- 


In  the  Yukawa  sector  the  /3  functions  are 


JY 


u,<f,e  _  Y     ^    /__L_^(1)       , 


c?t 


lg^2'^uAe^  (157^2)2  ^u,d,e 


P,/  ^  p) 


(A.8) 


where  the  one-loop  contributions  are  given  by 


17 


^i^^  =  -(Yu^Yu  -  Yrfty^)  +  F2(5)  -  (-^?  +  -A  +  S^i)  , 


20' 


j(l) 


/^y  =  2(Y/Yrf  -  Y„^Y.)  +  1-2(5)  -  (-9?  +  4^2'  +  %3)  , 
/3(^)  =  ^YetYe+F2(5)-^(5?+^i), 


(A.9) 


with 


Yl{S)  =  Tr{3Y«^Y«  +  SY^tY^i  +  Ye^Ye}  , 


(A.IO) 
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and  the  two-loop  contributions  are  given  by 

^?  =  |(YutY„)2  -  Y«tY,Y,tY,  -  iY,tYrfY„tY„  +  j(Y,tYj)2 

+  Y2{S){\Yd^Yd  -  -YJYu)  -  X4(5)  +  \\"  -  2A(3Y,tY„  +  Y^ty^) 

+  (^.?  +  f  .5  +  16.I)Y«tY.  -  (|,?  -  ^,i  +  16,i)Y,tY, 

5  9        29         4       9ool9o2.35  s4o22 

+  2^'^^'^^  +  ^  200  +  45"^^^'  "  20^^^2  +  Y^  AI  -  ( Y  -  "y)^2  +  9^253 
"404      80  '     . 

/3f  =  ^(YrftYrf)2  -  Y,tYrfY„tY«  -  ^Y.tY„Y,tY,  +  ^(Y„tY„)2 

+  Y2{S){\Yu'^Yu  -  ^YrftYrf)  -  X4(5)  +  ^a2  -  2A(3YdtYrf  +  Y,t Y„) 

5..  .^^         29         1  .       27   2  2      31    02        35  .4 

+  2^4(5)  -  (^  +  45"3^^1  ~  20^^^2  +  Y^5i53  -  (j  -  "3)fi'2 

"22      ,404       80  . 

4^^  =  |(YetYe)2  -  ^F2(5)YetYe  -  X4(5)  +  ^A^  -  GAYe^Ye 

'387   2      135   o^^ri^r       St^.^n      fsi        11      ,4      27   2  2 


with 


,35  4 


^^,2^  9    2^3   •2^^.   rv  t^ 


n(5)  =  (^?l  +  ^^2  +  8ff|)Tr{Y„TY,} 

+  (^^1  +  ^^2  +  8g|)Tr{YrftYrf}  +  ^(^2  ^  ^2)T,{YetYe}  , 


(A.ll) 


(A.12) 


(A.13) 


and 

X4(5)  =  ^Tr{3(Y«tYu)2  +  3(Y/Yrf)2  +  (YetYe)^ 

-  ^Y„tY,YrftY^}  . 
In  the  Higgs  sector,  the  /3  functions  for  the  quartic  couphng  and  the  vacuum 
expectation  of  the  scaleir  field  axe  presented.  Here  a  discrepancy  in  the  one-loop 
contribution  to  the  quartic  coupling  of  Ref.  7  is  corrected 

^  =  ^tf  ^  + -Vl/^f  '  (A.14) 
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where  the  one-loop  contribution  is  given  by 

+  4F2(5)A-4if(S), 
with 

H{S)  =  Tr{3(Y«tYu)2  +  Z{Yd^Ydf  +  (Ye^Ye)^}  ,  (A.16) 

and  the  two-loop  contribution  is  given  by  (discrepancies  found  by  Ford  et  al. 
[73]  in  Machacek  and  Vaughn    [7]  are  corrected), 

/?f  =  -78A3  +  18(^,?  +  3,|)A2  -  [  (^  -  10n,),|  -  ^^1 

9  ,229       50         .  .,        497      ^         «      3,97"    8         o  . 

"9  ,239       40         .2       27    59      40         a    ^ 

-  25^^  +  y"^)5l52  -  Y^(^  +  -^n,)g, 

-64^|Tr{(Yu^Y„)2  +  (Yrf^Yrf)2} 

-  g5r?Tr{2( Y„tY«)2  -  ( YrftY^)2  +  3( Ye^Ye)^} 

-  ^52^(5)  +  10A[  (^5?  +  ^52  +  8g|)Tr{Y„tY,}  (A.17) 
+  (^51  +  ^?2  +  8gi)Tr{YrftYrf}  +  ^(g2  +  ^2)Tr{YetYe}  ] 

+  \9l[  (-^<7?  +  2l5|)Tr{Y,tYu}  +  i\g\  +  9gi)Tr{Y/Y^} 
+  (-y^?  +  ll^|)Tr{YetYe}  ]  -  24A2f2(5)  -  \H{S) 
+  6ATr{Y„^Y„YrftY^}  +  20Tr{3(YutY„)3  +  3(Y/Yrf)3 
+  ( Ye^Ye)^}  -  12Tr{YutY,(Y„tY„  +  YdtY^)YrftY^}  . 
The  /3  function  for  the  vacuum  expectation  value  of  the  scalar  field  is 

dt     ~167r2^       +(167r2)2^       '  <^-^^^ 

where  the  one-loop  contribution  is  given  by 

i^'^  =  \{lo\  +  gl)-Y2{S).  (A.19) 
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and  the  two-loop  contribution  is  given  by 

^93      "l      ,4        511       5      ,4       27   2   2  ^^'^^^ 

These  expressions  were  arrived  at  using  the  general  formulas  provided  in  Ref.  7 


for  the  anomalous  dimension  of  the  scalar  field,  choosing  the  Landau  gauge. 

In  the  low  energy  regime  the  effective  theory  is  SU{3)c  x  t^(l)£'M-  The 
general  formula  of  Ref.  74  is  used  to  arrive  at  the  fS  functions  for  the  respective 
gauge  couplings 

+  15""  ^"'''S  '^•2'* 

,5033,  .325,  2      2857,    gl 


and 


de       .16  4  4     ,     e^  ,64  4  ,     e^ 

=  [-r-nu  +  -nj  +  -n/J— -^  +  [— n«  +  — n^  +  4n/]- 


c^f       ^9     "  '   9    '^  '   3    '^(47r)2   '  ^27    ""27    ■'^      '^(47r)4 
.64  16      ,  e^g^ 


(A.22) 


9     "^        9     '^^(47r)4  ' 
where  n^,  Hjf,  and  n/  are  the  number  of  up-type  quarks,  down-type  quarks,  and 

leptons,  respectively.  In  Eq.  (A.21),  the  three-loop  pure  QCD  contribution  to 

the  l3  function  of  g^  have  also  been  included    [75]. 

For  the  evolution  of  the  fermion  masses,  Ref.  76  is  used.  It  is  known  that 

there  is  an  error  in  their  printed  formula    [77].    This  typographical  error  is 

found  and  corrected.  The  calculation  of  the  two-loop  contribution  to  the  mass 

anomalous  dimension  in  QCD  using  the  corrected  formula  agrees  with  the 

result  obtained  in  Ref.  77.  Equation  (2.26)  of  Ref.  76  should  read 

3  32     • 

(A.23) 


76 

Using  the  corrected  expression,  the  following  mass  anomalous  dimension  is 
computed.  The  fermion  masses  in  the  low  energy  theory  then  evolve  as  follows 

dm 


dt 


=  l{l,q)^  '  (A-24) 


where  the  /  and  q  refer  to  a  particular  lepton  or  quark,  and  where 

m'?)=T'(/,<?)(i;^+^(/,9)(i;;]2 

I   .333    g3 

+  ^(9)  (47r)6  • 

The  superscripts  1  and  3  refer  to  the  ^/(l)^;^   and  SU{Z)(j  contributions, 
respectively.  Explicitly,  the  above  coefficients  are  given  by 

in  =  0 
■<t,)  =  -8 

^13  _  ^33  _  n 

7(/)  -  7(/)  -  0 

11  oM  fSO  20  20     ,^2  ^^'^■"^^ 


w.i)  =  -3^(/,g)  +  ty "« +  y"rf  +  y  "']^(/,?) 


^13  _      ,^2 
33  404      40, 

QQQ       2  140  o  ''216 

7(,)   =  gtyf  (""  +  '^df  +  (160C(3)  +  ~^){nu  +  rid)  -  3747]  , 

where  Qun)  is  the  electric  charge  of  a  given  lepton  or  quark,  and  (;(3)  = 
1.2020  ...  is  the  Riemann  zeta  function  evaluated  at  three.  In  the  mass  anoma- 
lous dimension  for  the  quarks  above,  the  three-loop  pure  QCD  contribution 
,333    r7c 


7,  V     [75]  have  also  been  included. 


APPENDIX  B 
CALCULATING  THE  MSSM  /3  FUNCTIONS 

In  this  appendix,  some  useful  results  in  calculating  the  j3  function  of  the 
MSSM  axe  included.  Also  some  sample  calculations  are  presented.  Although, 
to  one  loop  the  MSSM  0  functions  have  appeared  in  the  literature,  one  is  not 
aware  of  the  two  loop  ones'  appearance  except  in  some  approximate  form,  such 
as  keeping  only  the  contribution  of  the  heaviest  family. 

There  are  at  least  two  ways  to  proceed.  References  [78,79]  give  general 
formulas  valid  to  two  loops  to  compute  the  fS  functions  of  gauge  and  Yukawa 
couplings  in  a  supersymmetric  gauge  theory.  These  were  used  to  calculate 
the  two  loop  gauge  and  Yukawa  coupling  /?  functions  appearing  in  the  next 
appendix.  Reference  [7]  gives  formulas  valid  to  two  loops  to  compute  the  /? 
functions  of  gauge,  Yukawa,  and  scalar  quartic  couplings  in  a  general  gauge 
theory.  For  the  purpose  of  including  thresholds,  this  approach  is  more  useful, 
and  it  is  this  manner  that  the  following  examples  are  done.  The  results  of  this 
second  approach  were  checked  against  the  first  one,  and  they  agreed  (as  they 
must). 

Not  all  calculations  are  included,  as  this  would  be  useless  and  not  enlight- 
ening. However,  it  is  hoped  that  the  examples  will  give  some  flavor  of  the 
endeavor  and  that  the  results  will  prove  useful  to  anyone  wishing  to  pursue 
such  calculations  further. 
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B.l  Gauge  Couplings 

B.1.1  Group  Invariants 
Tables  6  and  7  will  be  useful  in  the  calculations  of  the  ^  functions  of  the 
MSSM.  Also  useful  are  the  following  group  invariants 

C2{R)  =  T^(i2)T^(iZ)  =  C2iR)  1{R)  .  (B.1.1.1) 

Using  this  definition  and 

Tr{T^(R)T^(R)}  =  T2(R)<5^^  (B.l. 1.2) 

the  following  identity  is  derived 

Tr{C2(R)}  =  diG)T2{R) 

C2iR)Tv{l{R)}  =  rT2{R)  (B.1.1.3) 

C2{R)d{R)  =  T2{R)r  , 
where  r  =  d{G)  is  the  rank  of  the  group,  and  8       =  d{G).  From  this,  it  follows 
that 

C2{G)  =  T2{G)  .  (B.l. 1.4) 

The  Yukawa  interactions  arising  in  the  SUSY  Lagrangian  can  be  written 
4^uklr.a  =  -^^9[<l>*i  T,j  ^l^jX^  -  i'j  Tfi  <i>,^^]  ,  (B.1.1.5) 

where  T  are  the  group  generators,  and  overlined  fermi  (Grassmann)  fields 
are  right-handed.  Now  make  the  correspondence 

„     SUSY      -  ^       A 

^  ^^'^  ^  (B.l. 1.6) 

Y -^(0  =  +^-^^^1  • 

Written  this  way,  the  i  (i.e.,  (t>[i))  subscript  represents  the  i     group  component 

of  the  scalar,  and  the  j  represents  the  group  component  of  the  fermion.  The  i 
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encodes  the  scalar  information  in  SUSY  whereas  the  a  does  so  in  the  notation 
of  Ref.  7.  When  taking  a  Yukawa  term  from  Ref.  7  and  "supersymmetrizing"  it, 
the  trace  need  not  represent  a  trace  unless  aJl  scalar  indices  (i.e.,  the  superscript 
on  the  Y's)  are  contracted.  For  example,  in  the  one- loop  contribution  to  the 
scalar  anomalous  dimension 

TrlY^Y^t}  =  [-iV2gT(l][iy/2gT^.] 


=  2fir  C2ij 
=  •2g^C2{S)6ij 
Now  to  compute  Y^(F). 

Tr{C2(F)Y^Y^^}  =  2g'^C2{F)TT{  I  T^T^} 
=  2g'^C2{F)TT{T^T'^]S^^ 
=  2g^C2{F)T2iR)6^^S^^ 
=  2g^d(G)C2iF)T2{R) 


(B.1.1.7) 


(B.1.1.8) 


y^iF)  =  -^Tt{C2(F)Y^Y^^} 

diG)  (B. 1.1.9) 

=  2g^C2{F)T2{R)  . 
These  results  may  be  used  to  derive  the  SUSY  gauge  beta  function.    In  the 

general  (but  single  simple  group  G)  case,  at  one  loop 

(47r)V^V5^  =  It2{F)  +  ^T2{S)  -  jC2{G)  ,  (B.1.1.10) 

where  F,  5,  and  G  stand  for  the  fermion,  scalar,  and  adjoint  representations, 
respectively.  To  go  to  the  SUSY  from,  the  following  correspondences  are  used 

F yR  +  G  . 

S yR  . 
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Whence 

(4^)^41/5y /^^  =  li^^iR)  +  C2(G)]  +  It2{R)  -  yC2(G)  ^  ^  ^ 

=  T^iR)  -  3C2(G)  . 
At  two  loops,  the  only  extra  complication  is  a  Yukawa  contribution  in  the 
general  case  that  must  be  supersymmetrized  as  discussed  above 

(47r)4/?(2)//  =  [jC2{G)  +  2C2{F)]T2{F)  +  [^C2(G)  +  4C2(5)]T(5) 

'?4 

-  jC2{Gf  -  Y^{F)  . 

(B. 1.1. 13) 
Applying  the  stated  correspondences  yields 

(47r)2/3g)^.//  =  jC2{G)[T2{R)  +  T2{G)]  +  ^C2(G)T2(i?) 

+  2[C2iR)T2iR)  +  C2{G)T2{G)]  +  AC2{R)T2{R) 

-  ^C2{Gf  -  [Y^{R)  4-  F4(G)] 

=  4C2(G)T2(i?)  +  QC2{R)T2{R)  -  6C2(G)2 

-[2C2{R)T2iR)  +  2C2{G)T2{R)] 

=  2C2iG)T2iR)  +  4C2{R)T2{R)  -  6C2{Gf  . 

(B. 1.1. 14) 


B.1.2  Fierzing 
A  matrix  M  can  be  expressed  in  terms  of  the  group  G's  generators 

M^f-'^T^  +  fl.  (B. 1.2.1) 

Taking  the  trace  of  both  sides  yields 

Tr{M}  =  f^d(R)  .  (B. 1.2.2) 

Multiplying  by  T°  and  taking  the  trace  gives 

Tr{MT^}  =  f^T2(R)<5'^^  .  (B. 1.2.3) 
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Hence 

M  =  -I^TrlMT^jT^  +  :J7^Tr{M}  I  .  (B.1.2.4) 

T2(R)  d(R) 

If  M  =  $$t^  it  follows  that 

**^  "  T^)^^^^^^^^''^  +  d[R)^^^^^  ^  '  (B.1.2.5) 

where  use  has  been  made  of  Tr{$$tT^}   =  Trj^^T^^}   =  $tT^$.     Now 
multiplying  on  the  left  by  $^  and  on  the  right  by  $  gives 

i^^^f  =  -^($tT^$)($tT>l$)  +  J_($t$)2  (B.1.2.6) 

Combining  like  terms  yields 

($tT^$)($tT^$)  =  T2{R)[l  -  3;^]($^$)2  .  (B.1.2.7) 

d(K) 

In  the  fundamental  representation  of  SU{n),  T2{R)   =   1/2  and  d{R)   =  n. 
Taking  SU{2)  as  an  example,  with  $  a  doublet,  it  follows  that 

($tx^$)($tx^$)  =  i($t$)2  _  (B. 1.2.8) 

It  is  also  useful  to  have  the  analogous  expression  for  the  case  when  there  are 
two  different  SU{2)  fields  involved 

(B.1.2.9) 
which  follows  immediately  from  Eq.  (B.1.2.5). 
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B.2  Anomalous  Dimension  of  the  Scalar  Field 
Using  the  Machacek- Vaughn  convention,  the  Yukawa  sector  of  the  SUSY 
Standard  Model  is  written 

ciE^^Ll  +  ^L^^dLi  +  tiE^dLL  +  c.c. 

=  ULia'^ij^uAii(^2)ABQLjpB^al3  +  KA(^^2)ABQL*pB^ji^L*ia  +  " " "    • 

(B.2.1) 

The  zT2's  involved  in  the  Weyl  spinor  product  have  been  suppressed  and  so 
have  the  icr2's  involved  in  the  SU{2)  product,  although  the  iag's  are  displayed 
in  the  second  step.  The  SU{3)  product  follows  since  the  barred  fields  are  anti- 
quarks.  All  fermions  are  represented  by  Grassmann  fields.  The  c.c.  is  explicitly 

calculated  in  the  example  below 

[XLiir2)r]L]*  =  xl{iT2)*vl 

=  Xl{i'r2)ril 

=  [XR{i'r2)]{ir2}[{-iT2)r)i,] 

where  the  identity,  t/;^  =  ^^24'!^  has  been  used.  The  following  identifications 
can  be  made  from  Eq.( B.2.1) 


(B.2.2) 


yr  +  U{(f>A)  c  I   ■  \  TT 

^[.(.i>A)](Qjm  +^c,'M2)AB^ij 
^[{uia){<i>A)]»  ^  +^a0("^^>AB^ij 


Y 


-u((f>A)  ^      ..      .       ,^f 

{QjPB)[.iui^  -^pa{i<^2)BA^ji  , 

v^-u(uta)  .       ,.  ,,t 


^[idia).]iQjpE)+^c.0{i^2)AB^ij 
^it(0i)](Qi/?F) +<^«^( '^2  ).45  D,-_,- 

^[{dia)(<pA)]>  "  +^c.0{i<^2)AB^ij 


-u(Qj0B) 
>[{<t>A){um)] 


r-d(<i>A) 


=  -^Paii<^2)BA^j 


(B.2.3) 


{Qjl3B)Hdi^  -*/?a(^^2)fl^D^.. 

^-d{dia)  -      ,.       ,        „t 

^Qj/^5)[(M)iI-'^/?'>(^^2)fl.4Djf 

^-d{QjpB)     _      .        .  ^t 


(B.2.4) 
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,  ,+(/(ei)  ,  ,  .       ,        ^  Vr-(i(ez)  /  •       \       T7>t 

^[.{M)]{LjB)=  +(i<^2}AB^ij  ,  ^iLjB)[(M)*r  '^""^'BA^ji  , 

(B.2.5) 
The  +  and  —  on  the  Y's  denote  whether  or  not  the  $  is  complex  conju- 
gated, respectively.  The  upper  index  always  denotes  the  scalar,  and  a  bullet  is 
placed  in  its  unoccupied  lower  position.  Since  two  of  the  three  fields  are  either 
transposed  or  not  (in  the  spinor,  SU(2)  doublet.  SU{Z)  triplet,  or  generation 
sense),  they  appear  grouped  in  square  brackets.  Dots  over  Y's  indicate  that 
the  Yukawa  coupling  matrix  it  represents  is  daggerred. 
Consider  the  following  example  from  Ref.  7 

X  =  ^TrlY^'Y^^Y^Y^"}  +  Tr{Y^Yt<^Y^Yt<=}  ,  (B.2.6) 

a  quantity  that  enters  the  scalar  2-loop  anomalous  dimension.  In  the  above 
notation  (suppressing  some  indices  at  first),  the  first  term  reads  (in  the  $u 
case) 

^^i+A){-B)  ^  ^^^y+u(M)y-u('^B)yy}  ,  (B.2.7) 

but  to  this  must  be  added  another  term 

_^.^+.4)(-fl)  ^  ^^^y+u(M)y-"(<^B)yY}  ,  (B.2.S) 

since  there  must  be  alternating  daggerred  Y's  in  these  expressions,  and  both 
equations  (B.2.7)  and  (B.2.8)  above  satisfy  this  (the  underlines  on  two  Y's 
indicate  that  their  scalar  indices  are  to  be  contracted).  Furthermore,  the  result 
must  be  symmetric,  in  this  case,  in  .4  and  B,  so  the  next  two  terms  must  also 

be  included 

y(-'-^)(+^)  ^  'Pi.{Y~"^'^'^^Y"'""^'^^)YY} 

(B.2.9) 
■^{-A){+B)  ^  ^^,y-u(M)y+"('^B)yy| 
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Putting  it  ail  togetlier,  tlie  first  term  in  Eq.  (B.2.6)  is  given  by 

X  =  ^-[{X+-  +  X+-)  +  (Xf  +  +  X^+)]  .  (B.2.10) 

Similarly  the  second  term  in  Eq.  (B.2.6)  is 

y  =  Im  +  5^2)+- + (Yi + 12)-+] .  (B.2.11) 

Now  to  compute  these  contributions  (the  flavor  and  color  indices  are  suppressed 
since  these  contract  in  a  straightforward  way) 

-^1       -  -^^t  »u.(QC)  ^(QC)»u^'u.(QD) '(QD).u  ^  *u(<^D).    ^.(<^D)u    J-" 

=  3Tr{(U^U)2)(icr2)AC(-i^2)CB[  (i'72)ED(-i<72)DE  +  (if^2)DE(-i<72)ED  ] 
=  12Tr{(UtU)2}<5AB  , 

X}-  =  0  , 
Xf  +  =  0  , 

V-+  _  rr.   f^-u(<i>A)^+ui<pB)r  v^ -u(«iE) ^+u( <^E )       ^-di4>E)^+di4>E) 
2       -  ^^\  ^  (QC).u  ^  u.(QD)  ^   ^  (QD).u  ^  u.(QC)  ^  *  (QD).d  *  d.(QC) 

^-U{u)  -yr+uiu)  ^-d(d)  ^+<i(d)  1-, 

"^  "^  (QD)(<t>E)*     •(<i>E)(QC)  ^  *  (QD){4>E)*  *  .(«i^)(QC)  J-f 

=  2-3Tr{(UtU)2utUDtD}(-ia2)cA(i<72)BD(-i'^2)DE(i^2)EC 

=  6Tr{(U^U)2  +  utUD^D}(!)AB  • 

(B.2.12) 

Some  of  the  3's  come  from  the  trace  over  the  suppressed  color  delta  functions. 

Finally, 

X  =  3Tr{3(UtU)2  +  utUD^D}<5AB  •  (B.2.13) 


Similarly, 


5-;+-  _  ']r'i-|Y+"*'*^)YY~"^'^^^Yl  =  0 
Y~+  =  Ti-;y~"^''^^^YY''""''^^^Y}  =  0 
Y — ^  =  Xr{Y~"*'^'^^YY"^"*"^^^Yl  =  0  . 


(B.2.14) 
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The  full  answer  is  then 

Xu  =  ^X+Y=  ^Tr{3(UtU)2  +  utUDto}  .  (B.2.15) 


B.3  Yukawa  Couplings 
The  Yukawa  interactions  arising  from  the  D-terms  in  the  SUSY  Lagrangian 
can  be  written 

4ufLa  =  -^^9[<f>*i  Tif  rPjX^  -  ^j  Tjf  4>i'>^^]  ,  (B.3.1) 

where  the  overhned  fermi  (Grassmann)  fields  are  right-handed.  Now  make  the 
correspondence 

The  n's  a's,  ©'s,  etc.  represent  field-types  (i.e.,  Q,  u,  J,  etc.).  Note  that  the 
Y's  which  represent  T^'s  are  diagonal  in  type  (i.e.,  both  indices  are  n's  as 
exemplified  above).  The  -\-  and  -  indicate  that  the  scalar  index  is  complex 
conjugated  (equivalently,  is  the  left  subscript  on  Y's)  or  not  complex  conju- 
gated (equivalently,  is  the  right  subscript  on  Y's),  respectively.  Written  this 
way,  the  i  above  represents  the  i^^  group  component  of  the  scalar,  and  the 
j  represents  the  group  component  of  the  fermion.  Contracting  two  Y's  (i.e., 
using  underscores)  implies  in  their  case  contracting  the  superscript,  group  gen- 
erator indices  (the  i?'s). 


(B.3.2) 


(B.3.3) 


86 
In  order  to  construct  the  one  loop  SUSY  beta  function  formula  for  the 
Yukawa  matrices  the  following  must  be  supersymmetrized 

\  J[d*]a 

_~-iXR)~+(\R)    +{<i>)    ,  ^+{m^-i>^R)^+(<t>) 

Y^Y^Y    ^^^  JY+(<i>)Y^i>^R)Y'^i^R)] 
_^M4>)  ^-(XR)^+{^R) 

=Yi,(^^  C2(a) 

[D«]a     ^^    ' 

-^  2g"^Y('^)C2(i?)  . 
Furthermore,  in  obtaining  the  formulas  for  the  beta  functions  of  the  different 
Yukawas  the  "external  subscripts"  are  generally  two  fermions  (not  scalarinos) 
which  are  represent  by  n's,  a's,  q's,  etc.,  and  the  "external  superscript"  is  a 
scalar  (a  Higgs  not  an  sfermion)  which  is  represented  \iy  ^{<f)-  Hence,  it  should 
be  clear  from  inspection  that  the  following  terms  cannot  be  constructed  for  any 
choice  of  "internal  subscripts" 

YTrlY^Y^} >(Y±(^'^)Tr{Y^(^^^Y('^H)DA  ^  0  , 

~  (B.3.4) 

YY°ty    ^^^^  ^^Y^[\K)Y^<i>-)Y^{\R')^^^  ^  Q  _ 

Application  of  Eq.  (B.3.3)  to  the  up- Yukawa  (Yr    \iqd\\  see  section  where 
the  SUSY  Y's  are  defined)  yields 

=  (c2(u)-fC2($t,))Y^-;;;](g5)' 
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where  use  was  made  of  the  diagonal  nature  of  C2n^  ~  ^nAC2(n).   Note  that 
C2  on  the  left  can  act  on  both  left  subscripts.  Likewise 


(yC2{R)) 


+<<!>) 


[u.]{QB) 


(B.3.6) 


(B.3.7) 


(B.3.8) 


The  standard  one  loop  beta  function  for  the  Yukawas  is 
(47r)2/3''  =  ^  [yY^Y°  +  V'Y^y]  +  2YY''^ Y  +  YTr{Y^ Y*} 

Supersymmetrizing  according  to  the  results  derived  above  gives 
(47r)2/3°f;^  =  i  [yY^Y"  +  Y^'Y^y]  +  2YY"^Y  +  YTrjY^Y^} 
-{3-~6g)g^[c2{R).Y^}  , 
where  6g  equals  1  or  0  depending  on  whether  one  is  above  or  below,  respectively, 
the  mass  threshold  of  the  gaugino  in  question.  The  labels  Tj,  i  =  1, . . . ,  5,  will 
be  used  to  denote  the  five  terms  appearing  in  Eq.  (B.3.8)  when  computing  the 
supersymmetric  beta  function  for  the  up- Yukawa 
T^  =  (^YY^Y"*""^'^^') 


r  +  u{<i>A) 


[u»]{QC)     iQC)[*u]  ^  "  [ui4>C)]*  "  •[i4C)u]\  '■  [u.]{QB) 


(B.3.9) 


=  («<^2)^bU 


4Utu 


T2  =  ('y+"*'^'y^y) 


]{QB) 


=  Y 


-¥u(<i)A) 


u.\{QC) 
+  Y 


y.-u{<pD)  -y+u{4>D)         -y— u(u)  -^- 

[QC)[*u\     [u.](QB)  ^      (QOUD).]  ^  [.(4,D)]{QB) 


r  +  u(u) 


d{<i>D)  ^+d(<i>D)     ,   ^-d{d) 


(QC)[*d]     [d*](QB)  +  ^  (gC)[(^D).]  ^  [.i4>D)]{QB) 


;b.3.io) 


=  ('^2)/15U 


2U'U  +  2D^D 


[u.](QB) 
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=  0, 

^       V-     ^-  V[u.](QB) 

-  ^  [u.](QB)  ^^^  ^  (QD)[.u]  ^  [u.](QD)/ 
=  (^>2)C5UTr{3Utu}(-i<72)DC(i^2)AD 
=  (^>2)^5U[Tr{3UtU}]  , 

3 


k=l 
4        1         1      2   ,   ,n   ,    3    ,    3,   o 


=  (^-2).45U  [(g  +  ^  +  1)9'  +  (0  +  i  +  i)52 


4       4  9 


=  («<^2).4bU 


13   2       "2       ^2 
Y^9    +  ^92  +  353 


(B.3.11) 


(B.3.12) 


(B.3.13) 


In  T4  the  3  comes  from  the  trace  over  the  suppressed  color  indices.  Also,  in 
T5  use  was  made  of  the  fact  that  if  the  gauge  group  is  not  simple  but  a  direct 
product  Gi  X  •  •  •  X  Gn  with  couplings  51, . .  ■  ,?n,  then  g  C2  —>■  ^|t=i  91^2   ■ 


B.4  Thresholds 
To  implement  the  super  particle  thresholds  in  the  minimal  low  energy  super 
gravity  model,  the  renormalization  group  /?  function  must  be  calculated  in  a 
form  that  has  not  appeared  in  the  literature.  In  the  following  example,  the  one 
loop  l3  function  of  g^  is  considered.  Starting  from  Eq.  (B. 1.1. 12) 

(47r)24^V?3  =  It2{Fz)  +  It2{Sz)  +  lc2{G-i)6g  -  jC2{Gz)  ,       (B.4.1) 

where  F^,  53  refer  to  the  fermion  and  scalar  representations,  and  G3  =  SU{3). 
Also 


1  .     ^  >  Mg    : 

0    .         H<    Mg       , 


(B.4.2) 
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where  Mg  is  the  mass  of  the  gluino.  When  dealing  with  a  direct  product  group, 
like  SU{3)  X  SU{2)  X  U(l), 

r2(i?3)  ^  T2{Rz)d(Ri)d(R2)  .  (B.4.3) 

From  the  definition  of  T2{R)  in  Eq.  (B.1.1.2),  one  obtains  the  following  result 
in  the  SU{3)  case 

T2{Rz)  =  2(i)iVg  +  (i)iV^  +  (i)iV^  ,  (B.4.4) 

where  Np  equals  the  number  of  generations  of  particle  p.  This  result  is  valid 
for  both  fermion  {R  =  F)  and  scalar  {R  =  S)  representations.  The  notation 
will  be  such  that  N  =  n  for  the  particle  and  iV  =  n  for  the  SUSY  partner. 
Equations  (B.4.1)  and  (B.4.4)  lead  to 

(47r)2/?(^Vff3  =  li^u+nd)  +  -hQ  +  -hjj+-h^+26g-n  .  (B.4.5) 

It  has  been  assumed  that  uq  =  [uu  +  n^)/2.  Also  the  fact  that  left-handed 
and  right-handed  quarks  of  a  given  flavor  have  the  same  mass  implies  n^  =  n^ 
and  nj  =  n^.  Note  that  this  reduces  to  the  right  standard  model  result  when 
h  =  Sg  =  0,  and  to  the  right  supersymmetric  result  [fi  >  MsuSY )  when  n  =  3, 
8g  =  \.  Similar  formulas  are  calculated  for  gi  and  52 

,2,3(1)/  3       2   17  5  5  1  1  4  1 

(47r)  ^1    Ig,  =  -{~nu  +  -n^  +  -Ue  +  -n.)  +  -hg  +  +-n^+  _^- 

1-1,1,  1 

9   n)     •}         ''2      1  1  11 

(47r)  /?^  7^^  =  -y  +  2("«  +  ""d)  +  g(ne  +  n;.)  +  -ng  +  -n^ 

1        "  1  4~- 

(B.4.6) 
For  the  gauge  couplings,  the  two  loop  contributions  were  also  calculated  in  this 
manner.  For  the  Yukawa  couplings,  this  form  of  the  /i  function  was  calculated 
to  one  loop. 
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Table  6:  Group  Theory  Factors 

T2{R) 

C2{R) 

diR) 

Ui 

SU2 

SU^ 

Ui 

SU2 

SU^ 

Ui 

SU2 

5C/3 

Q 

i^f 

1 
5 

1 
5 

i+kf 

3 

4 

4 
3 

2 

3 

u 

i-if 

0 

1 
5 

i-h' 

0 

4 
5 

1 

3 

1 

(+i)' 

0 

1 

5 

m)' 

0 

4 
5 

1 

3 

L 

i-i)' 

1 
5 

0 

(A)' 

3 

0 

2 

1 

e 

(+1)2 

0 

0 

(+1)2 

0 

0 

1 

1 

^u 

mi' 

1 
5 

0 

i+hj 

3 

0 

2 

1 

^d 

(A)' 

1 

7 

0 

(-f)' 

3 

4 

0 

2 

1 

Table  7 

G 

C2{G) 

diG) 

Ui 

0 

1 

SU2 

2 

3 

SUs 

3 

8 

APPENDIX  C 
THE  MSSM  l3  FUNCTIONS 

Using  some  of  the  notation  of  Falck    [80],  the  superpotential  and  soft  sym- 
metry breaking  potential  are  as  follows 

W  =  flYu^uQ  +  ^d^dQ  +  ^Ye$jL  +  fi^u^d  +  h.c.  , 
Vsoft  =  '^l^Uu  +  ^i^i^d  +  Bil(<^u^d  +  h.c.) 

+  mnQ' Q  +  rrij^L' L  +  m:^u  u  +  m^d  d  +  mgC  e 

+  Y,i  AiJvi'jli^uQj  +  A'JY'Jdi^iQj  +  j^Yi^^^dLj  +  h.c.  )  , 
hi 

1    ^ 

"  /=1 

(C.l) 
Various  (T2's  have  been  omitted,  and  a  sum  over  the  number  of  generations  is 
implied  in  the  squark  and  slepton  mass  terms.  Also,  hats  imply  superfields 
and  tildes  the  superpartners  of  the  given  fields. 

First  the  gauge  couplings 

dgi  _         1 


dt  167r2 

,3 


„3 


(C.2) 
where  t  =  In//  and  /  =  1,  2,  3,  corresponding  to  gauge  group  SU{'i)(j  x 
SU(2)i  X  U{\)y  of  the  Standard  Model.   The  various  coefficients  are  defined 


91 


to  be 
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■l-2n, 

0 


bl    —    —-    —   ^llg    , 

62  =  5  -  2ng  , 


63  =  9  -  2n 


9  ' 


/  38       6       88  \  /  Q 


(blk) 


I      14      8 

I        11  q  68        I 


Un    + 


55 
3 


and 


with  Ug  =  ^rifi. 


(2t)       14       i»  \ 
I       g       I        ,with/ 
4       4       0/ 


9  n 

-17       0 
0        -54 


=  K  ,  rf  ,  e 


:c.3] 


(C.4) 


(C.5) 


In  the  following,  the  beta  functions  for  the  parameters  of  the  superpotential 
are  listed. 


^  =  I^[  Tr{3Y„tY,  +  3Y^tY^  +  y^ty,}  -  3(^5?  +  g^) 


(C.6) 


In  the  Yukawa  sector  the  /3  functions  are 


dV 


u,d.e 


.,(1) 


,(2) 


dt  ^^,dJ-^Q^2^u,d,e+  ^iQ^2y2^u,d,e)  '  (^.7) 


where  the  one-loop  contributions  are  given  by 


13   o 


16 


,4^'  =  3YJYu  +  Yd^Yj  +  3Tr{YjY„}  -  (-gj  +  3gl  +  yy|)  , 

/?i^)  =  3Y/Yrf  +  Y,tY,  +  Tr{3YrftY^  +  Y^tY,}  -  {^gj  +  Zg\  +  ^^32)  , 


15^ 


9 


4^'  =  3YetYe  +  Tr{3Y/Yrf  +  Y^tYe}  -  (-5?  +  Zg\)  , 


(C.8) 
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and  the  two-loop  contributions  are  given  by 

;?i2)  ^  _4(Yu^Y,)2  -  2(Y/Y^)2  -  oY^^d'^jYu  -  9Tr{Y,tY.}Y«tY„ 

-  Tr{3Y/Y^  +  Ye^YelYrftY^  -  3Tr{3(Yu^Y,)2  +  Y^' YrfY«^Y„} 

+  {lgl  +  QglwJ^u  +  qghV^d  +  i^9l  +  16gi)Tr{Y«tY,} 
,26  403,   4      ,^  21,    4      ,32  304     4 

9  o       136   90      „  '?  1 

+  9i92  +  j^sfgi  +  8gM 

I3f  =  -AiYd^Yif  -  2{YjYuf  -  2YjY„YrftYd  -  3Tr{Yut Y„}Y^^Y„ 

-  3Tr{3Y/Yrf  +  Ye^YelY^tYrf 

-  3Tr{3(YrftYrf)2  +  (Ye^Ye)^  +  YrftY^Y^tY,}  +  ( ^gl)Yu^Yu 
+  (^51  +  6gi)Yd^Yrf  +  {-^l  +  16g|)Tr{Yd^Yrf} 

+  (^^?)Tr{YetYe}  +  (^n,  +  ^)g\  +  (6n,  -  y)^^ 

^  32  304     4        99      800      o22 

+  (y"3  -  -^)53  +  ?192  +  9^153  +  8^293 

/3f  ^  =  -4(Ye^Ye)2  -  3Tr{3YrftY^  +  Ye^Ye}Ye^Ye 

-  3Tr{3(Y^tY^)2  +  (Ye^Ye)^  +  Y^^Y^Yu'^Yu] 

+  (6g|)YJYe  +  (^g?)Tr{Y,tY}  +  (-^g?  +  16g|)Tr{Y/Yrf} 
,18  27^  4      ,^  21,  4      9   2  2 

(C.9) 
The  evolution  of  the  vacuum  expectation  values  of  the  Higgs's  is  given  by 

^^^^u,0d  _  J_    (1)        ,_L_J2) 

dt         ~  167r2^<*u,0d  ^  (I67r2)2^<?iu,0d  '  ^^-'^^ 

where  the  one-loop  contribution  is  given  by 

-'■i?  =  H^'  +  ^2)-3Tr{Y.tY,}, 


^i?  =  i^^^l  +  ^2)  -  3Tr{Y,tY,}  -  TrjYetYe} 


;c.i2) 
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and  the  two-loop  contribution  is  given  by 

ff^  =  ^Tr{3(YjY„.)2  +  3Y,tY„YrftYd} 

-(^^l  +  ^^2  +  20g|)Tr{Y„tY„} 

,279       1803         A      ,207      357         .        '27        9       .32 

7^^^^  =  ^Tr{3(Y/Yrf)2  +  3YrftY^Y«tY„  +  (YetYe)^} 

-  (f  5?  +  ^^1  +  20gi)Tr{Y,tY,}  -  (|g?  +  ^^i)Tr{Y,tYe} 

279   ,    1803      ^  4      ^207      357         4      ^7        9  22 

-^800+l600"^^^l-^"^  +  ^"^^^2-(^  +  ^n.)^l52- 

The  renormalization  group  equations  for  soft  symmetry  breaking  terms  are 

dt  IQir^'  Y^J  Ye^  1; 

+  2(aJ"'|F/-'"1^  +  3a5"^|f/"^|2)  -  6(^9fMi  +^|M2)  ]  , 

^  =  7^[  4(Y,Y,t)^-^4^^  +  5A;f -^(Y,tY,)^^-  -  3-4(Y,Y,tY,)^-^- 

^d  ^d  ^d 

+  [Af  -  a'J){Yu^Yu)''^%  +  2{YdYu^y'--Ai^% 

^d  ^d 

+  2(Aj-|Fe''"P  +  3a5-|F/-|2)  -  ll^J^i  -  6g|M2  -  y^lMg  ]  , 
^  =  T^[  4(Y,Y,t)=-^4S^\  +  5A|f  iiW(Y.tY„)^-i  -  3:%(Y„Y,tY,)'i 

+  (A?  -  Ai/)(YrftYrf)^-^'^  +  2(Y„Y^t):^^45;^ 


y'J  -  "    y 


(C.13) 
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(7777.  1 


dmi. 
M_ 

dt 


5 
dt 


+  IdlTviYm^}  -  ^gl\lf  ]  , 


drrij^ 


=  ^[  E  i^e^'i  VL  + -i  + -?.  +  i^m 


J 


A^^2rp^/v_2i       L2A.r2       o„2,r2 


J 


drriy, 
dt 


^^Q, 


|,?Tr{Fm2}  -  1|,2,,2  _  |^|^,2  ^ 


'J 


(C.14) 
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where,  as  in  Ref.  SO.  sums  are  implied  over  all  indices  not  appearing  on  the 
left  hand  side  and  where 

ng 

Tr{Fm2}  =  ^(m^_  -  2ml  +  ^|  "  ^l,  +  ^l)  ■  (CIS) 

i=l 

The  gaugino  masses  evolve  as  follows 

c/lnMi  1    ,     9 


APPENDIX  D 
THE  i3  FUNCTION  OF  THE  VEV 

Here  one  examines  under  what  conditions  the  /?  function  of  the  vev  agrees 
with  that  of  the  anomalous  dimension  of  the  corresponding  scalar  filed  by 
considering  a  simple  toy  model.  Spontaneously  broken  scalar  electrodynamics 
is  considered.  The  following  bare  Lagrangian  describes  the  theory 

where  D^  —  d^,  -  ieoQAg^  is  the  covariant  derivative.  $o  =  {(f>()  +  ixo)/V2  is 
a  complex  scalar  field  of  charge  one  (Q^g  =  *o)-  The  form  of  its  potential  is 

F($o)  =  Aol^J^o  -  Ifif  ,  (D.2) 

where  /q  is  the  minimum  of  the  bare  potential.  If  vq  represents  the  bare 
vacuum  expectation  i-alue  of  the  scalar  field,  $o»  then  to  expand  around  the 
true  vacuum  the  identification  is  made,  $o  =  "^o  +  ^-'O)  where  <  ^L  >=  0.  If 
the  primes  are  dropped,  this  transformation  amounts  to  taking 

<?>0  -^  "^O  +  VQ  , 

(D.3) 
xo  -*  \0  • 

In  terms  of  these,  the  gauge  fixing  (GF)  and  ghost  Lagrangians  are  given  by 

^  1    _  2 

^GFO  =  -T^id-  Aq  +  60^0^0 Yo) 

1  .  1   .    .  (D-4) 

=  -^(^  ■  --^O)^  +  eoVQAo  ■  d\o  -  -elvl^oxo 


and 


^ghosta  =  ^o(7-)^0  .  (D.5) 
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where  F  =  d  ■  Aq  +  eoyoifoXO  is  the  gauge  constraint,  a  is  the  gauge  function, 
and  770  is  the  ghost  field.  Under  a  gauge  transformation 

6(i>Q  =  CQaxQ  ,  (D.6) 

^XO  =  -eoOii(l>o  +  vq)  . 

Therefore 

6F  =  -Da  -  eQVQ^Q{(l)Q  +  VQ)a  , 

SF  ,  ^°-^^ 

■7-  =  -n  -  eoUo^o(</'0  +  ^^o)  • 

da 
Consequently  the  ghost  Lagrangian  is 

^0(^)^?0  =  d^iVod'^VO  -  ^oeo"0'?0^0  "  eot'o^O'70^0'?'0  (D.8) 


6 


a 


In  terms  of  (/iq  and  xo  the  covariant  derivative  term  above  can  be  rewritten 

l^<^0l^  =  Id^hd^h  +  \d^,xod^Xo  +  eoAo^ixod''<t>o  -  <f>od^'xo) 

+  ^el^Uo  ■  Ao  +  \elxlAo  '  -^o  (D.9) 

J.     o     9  o 

+  ^^O^O^O  •  ^0  +  eoVo(l>oAQ  ■  Aq  -  cqvqAq  ■  dxQ 
The  potential  becomes 

V  =  -Ao0o  +  ^o''o4  +  :^^o(3j;o  -  fh^o  +  ^oM4  -  fo)^ 


^  -       ,4         Ix     /„2         r2x,  2 


+  ^kxo  +  7,ki4  -  fo)\o  (D.IO) 

loo  0 

The  following  are  some  cf  the  relations  amongst  the  bare  and  renormalized 

parameters  that  will  be  needed  later 

1 

(D.ll) 

h  =  zjf, 


\o-- 

2* 

VQ   = 

i 
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with 

Z<^  =  1  +  8Z^  , 

Za  =  1  +  bZx  , 

(D.12) 

Zv    =    I   +  SZy     , 

Z/  =  1  +  (^'Z_/:  . 
The  6Z  terms  axe  considered  small  and  all  calculations  are  done  to  first  order  in 
these.  Furthermore,  the  following  relations  between  the  bare  and  renormalized 
n-point  functions  will  be  useful  later 

rSi'^  =  4[-V'o(^o-/o)  +  To], 

r^2)  =  z^[  p2  -  Ao(3t;o^  -  fi)  +  So(p2)  ]  , 

3  (D.13) 

r^^U  z|  [-6Ao^o  +  eo], 

rji'^  =  4^  -6Ao  +  Ao  ]  , 

where  Tq,  SqCp  )5  ©Oi  ^^^^  Aq  are  the  connected  loop  contributions  to  the  one, 
two,  three,  and  four  point  functions,  respectively.  From  above  the  following 
identifications  are  made 

Ml  =  elvl  , 

ml  =  Xo{3vl  -  f'l)  , 

(D.14) 
rnxl  =  ^Oivl  -  fS)  +  ^oM'i  ' 

mr]Q  =  ^qMq    , 
where  Mq,  mo,  ^t^q,  and  mrjQ  are  the  masses  of  the  gauge  boson,  Higgs  boson, 
Nambu-Goldstone  boson,  and  ghost,  respectively.    Making  the  choice  v  =  / 
(Note:  This  is  true  at  one  renormalization  mass  scale  only)  and  dropping  the 
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subscript,  the  first  of  Eq.  (D.13)  becomes 


1 


'^\2v 


r^^^  =  ^|[-A^'^(t''^.-/^^/)  +  ro], 


z 


1 


^\  '  '  ■  (D.15) 


=  Zl[-\v^-^{6Zy-6Zf)  +  TQ 


1 


=  _A„3:?i|i(«Z„  _  SZ,)  +  T  , 

=  -Xv^(6Zv  -  6Zf)  +  T  . 
This  leads  to  the  result 


T 

6Zy  -  6Zf  =  -^  .  (D.16) 

Similarly  for  the  second  of  Eq.  (D.13)  it  follows  that 
r(2)  =  Z4  p2  _  A||(3t;2z,  -  f^Zf)  +  Vo(p2)  ]  , 

=  Z^[  p2  _  2At;2||(l  +  |<5Z,  -  IsZf)  +  i:o(p2)  ]  , 
=  Z^p2  _  ^2|A(i  +  3^^^  _  1^^^)  ^  3,(^2)  ^ 

(D.17) 
Rewriting  S(p2)  =  v;(0)  +  !]'(0)p2,  where  the  prime  denotes  differentiation 
with  respect  to  p^,  Eq.  (D.17)  implies 

6Z^  =  -^'{0)  (D.18) 

and 


Q  -1 

=  p2  _  ^2  ^  ^^^^2  _  ^2^^2x  -  6Z^  +  -SZy  -  -SZf)  +  E(p^)  . 


3  1  ^(0) 

6Zx  -  SZ^  +  -8Zy  -  -SZf  =  -^  (D.19) 


Above  the  identification  has  been  made 


m'^  =  2Au-  (D.20) 


101 

for  the  renormalized  mass  of  the  Higgs  boson.  Now  for  the  third  of  Eq.  (D.13) 

1 

r(3)  =  z|[-6A.^  +  0o], 

=  -6At>  -  6Ati(iZj  +  -<SZ„  -  -iZ^)  +  e  , 
From  Eq.  (D.21)  the  following  result  is  obtained 

6Zx  +  \{SZy-8Z^)  =  -^^.  (D.22) 

Finally  the  fourth  of  Eq.  (D.13)  gives 

r(4)  =  z2[-6A|  +  Ao], 

=  -6AZa  +  A,  (D-23) 

=  -6A  -  Q\8Zx  +  A  , 
from  which  follows  the  identification 

bZx  =  A  .  (D.24) 

The  symbols  T,  E(0),  Il'(O),  0,  and  A  represent  the  renormalized  loop  con- 
tributions to  the  various  n-point  functions.  They  represent  C(A  ,Ae  ,e  )  cor- 
rections to  the  tree  level  results.  They  are  correctly  calculated  in  terms  of 
renormalized  quantities  to  the  order  that  is  being  employed.  Consequently  the 
following  renormalized  masses  will  be  needed 


A/2 

2    2 

=  e  V    , 

m^ 

=  iM'^  . 

m^ 

=  ih'f  . 

(D.25) 
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The  various  bare  propagators  required  to  perform  the  calculations  are 


—  I 


(a)  the  gauge  boson  :      — — -[  g^^,  +  (^q  -  1),^    ^ /,   .  ] 

(b)  the  Higgs  boson  : 


k'^  —  m 


(c)  the  Nambu  —  Goldstone  boson 

— ? 

(d)  the  ghost 


(D.26) 


/t2 


'"XO 


k^  —  m 


2 


The  renormalized  propagators  are  arrived  at  by  simply  substituting  in  them 
the  proper  renormalized  masses. 


The  bare  couplings  are 

(i)  H-NG-7  : 
(ii)  2H-27  : 
(iii)  2NG-27: 
(iv)  H-27  : 
(v)  4H  : 
(vi)  3H  : 
(vii)  IH  : 
(viii)  4NG  : 
(ix)  2H-2NG  : 
(x)  H-2NG  : 
(xi)  H-2G  : 


,    [in)  (out). 

eo(Pi        -P2        ) 

-iXqvq  , 
-i\ovq[vI  -  f^)  . 

-iXQVQ  , 


(D.27) 


where  H,  NG,  7,  and  G  stand  for  Higgs,  Nambu- Goldstone,  photon,  and  ghost 
lines,  respectively. 
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The  results  of  the  calculations  are 

A 


!:(0)  ,       1     2      9e^,      1 

m^  2  2  \    loTT^e 

W(0)  =  -e2(3-O-V 

-^  =  (6A  +  ee^  +  3  y)t7:V  (°-28) 

Au'^  A    IbTT'^e 

--  =  (10A-ee2  +  3^)— ^ 
6Ai;  A    iDTT-'e 

A  _.       ,.  9        e^.      1 


^  =  (10A- 2^6^  +  3 


6A       '  ^^  A'l67r2e 

Using  Eqs.  (D.16),  (D.IS),  (D.19),  (D.22),  and  (D.24)  the  following  result  is 
obtained  for  the  counter  terms 


6Z^  =  e\z  -  0 


167r2f 


e\      1 


<5Z;,  =  (10A-2ee2  +  3^)— ^ 

^    ^^'^^  (D.29) 

0         e4        1 
^Z    =(_6A  +  3e^-3-)— ^ 

Since  the  /?  function  comes  from  these  pole  terms     [81],  it  is  evident  that 
j3v  —  i^ (f,  '"^^  ^  —  0,  the  Landau  gauge. 


APPENDIX  E 
NUMERICAL  TECHNIQUES 

The  Runge-Kutta  method  is  used  to  numerically  integrate  the  /?  functions. 
There  are  eighteen  coupled  first  order  differential  equations  involved  in  running 
the  Standard  Model  couplings.  At  one  loop  some  equations  decouple  from  the 
rest.  For  example,  the  gauge  couplings  are  decoupled  at  one  loop.  However, 
the  Yukawa  /?  functions  depend  on  both  gauge  couplings  and  Yukawas  even  at 
one  loop  as  in  Eqs.  (A. 8)  -(A.9).  One  is  working  with  two-loop  /?  functions  and 
these  are  all  coupled.  Standard  Runge-Kutta  programs  are  readily  available; 
however,  these  usually  assume  knowledge  of  initial  values  of  all  couplings  at 
the  same  scale  to,  where  t  will  denote  the  scaling  parameter.  This  presents 
somewhat  of  a  problem  since  different  couplings  may  be  known  at  significantly 
different  scales. 

The  method  employed  to  solve  the  initialization  problem  is  often  referred 
to  as  "shooting"'  [82].  Simply  put,  a  guess  is  made  for  the  initial  values  of 
all  the  parameters  at  the  chosen  scale  to-  The  parameters  are  then  evolved  to 
the  various  scales  where  one  has  known  values,  and  the  merits  of  the  guess  are 
assessed.  The  procedure  is  optimized,  and  one  thereby  arrives  at  a  solution. 
Arriving  at  initial  values  of  the  couplings  at  the  same  initial  scale  therefore 
involves  the  solution  of  a  possible  eighteen  nonlinear  equations  in  eighteen 
unknowns.  That  is.  in  the  most  general  case,  if  it  is  assumed  that  none  of 
parameters  are  known  at  the  desired  initial  scale  to,  then  eighteen  coupled 
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functions  of  eighteen  variables  can  be  defined  as  follows 

Fi(x(to))  =  RK{y:{to); if }i  -  Xi  ,  (E.l) 

where  'x.(to)  is  an  eighteen  component  vector  denoting  the  unknown  parameters 
at  to,  Xi  is  the  known  (possibly  experimentally  determined)  value  of  the  zth 
parameter  at  some  scale  tj,  and  RK{x{to);  ti}i  is  the  value  of  the  ith  parameter 
resulting  from  numerically  integrating  the  l3  evolution  equations  to  a  scale  tj 
given  initial  values,  x(to)  at  to-  Solution  routines  that  solve  N  simultaneous 
nonlinear  equations  in  N  unknowns,  that  is,  they  solve 

Fj(x«o))  =  0,  (E.2) 

are  then  used  to  find  x(to).  In  the  present  case  Xi{to),  i  —  1, ...,  18,  represent 
the  eighteen  parameters  of  the  Standard  Model,  and  the  initial  scale  i^o  = 
e  0  =  Mz  has  been  chosen.  The  task  has  been  appreciably  simplified  since  one 
has  most  initial  data  at  M^ .  Only  the  fermion  masses  are  taken  at  a  different 
scale.  This  reduces  the  number  of  simultaneous  equations  and  therefore  the 
computing  task. 

As  the  top  quark  and  Higgs  boson  masses  are  unknown  in  the  Standard 
Model  at  present,  in  the  process  of  the  analyses  one  is  free  to  choose  values 
for  these  masses  at  Mz  and  then  proceed  to  study  the  consequences.  In  other 
models  in  which  there  may  be  certain  constraints  [e.g.,  in  some  GUTs  the  b 
and  T  Yukawa  couplings  are  equal  at  the  scale  of  grand  unification),  one  may 
incorporate  these  constraints  into  the  functions  (E.2).  The  freedom  to  choose 
a  \-alue  for  an  unknown  parameter  may  be  replaced  by  such  a  constraint,  and 
this  may  result  in  a  definite  prediction  for  that  parameter.  Constraints  from 
grand  unification  and  supersymmetry  were  used  in  Ref.  3  to  arrive  at  possible 
values  for  the  top  quark  and  Higgs  boson  masses. 
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After  ail  data  is  obtained  at  Mz  by  employing  the  initialization  procedure 
described  above,  the  Runge-Kutta  routines  are  used  to  evolve  the  parameters 
to  any  mass  scale  /i.  Several  figures  are  provided  in  Section  3.7  displaying 
the  evolution  of  many  of  the  parameters  of  the  Standard  Model  for  the  case 
mt{Mz)  =  95  GeV.  As  discussed  in  Section  3.4,  given  the  evolution  of  the 
running  parameters  m<  and  A,  the  physical  top  quark  and  Higgs  boson  masses 
may  be  found  by  solving  Eqs.  (3.4.1)  and  (3.4.2),  respectively.  These  equations 
axe  solved  using  the  nonlinear  equation  solution  routines  described  above  in  the 
initialization  procedure.  For  example,  consider  the  above  situation  in  which 
mt{Mz)  =  95  GeV,  and  suppose  one  wishes  to  find  the  corresponding  physical 
mass  Mt  in  this  case.  The  solution  algorithm  may  be  described  as  follows:  A 
guess  is  made  for  Mt,  then  the  running  parameters  are  evolved  using  Runge- 
Kutta  to  this  mass  scale  {i.e.,  ^  —  Mt).  The  guess  is  tolerated  depending  on 
how  accurately  one  wishes  Eq.  (3.4.1)  to  be  satisfied  when  values  for  a.s{Mt) 
and  mt{Mt)  are  substituted.  The  solution  routines  effectively  optimize  this 
shooting  or  guessing  procedure  and  yield  a  value  for  Mt  corresponding  to  the 
value  mt{Mz).  In  this  particular  case,  the  result  is  Mt  =  100  GeV. 


APPENDIX  F 
EXPLICIT  FROM  OF  6{fi) 

In  Ref.  34  the  radiative  corrections  term  6{fi),  from  Eq.  (3.4.2),  is  derived. 
In  this  appendix,  its  expUcit  form  is  presented  as  it  appears  in  this  reference 
except  for  some  minor  notational  changes.  In  the  following,  s  and  c  refer  to 
sm9\y  and  cos^jy,  respectively.  Also,  i^  is  defined  to  be  the  ratio  MJj/M^. 

G    M^ 

<^(/^)  =  ^^u/i(e,/i)  +  /o(e,M)  +  rV-i(^/^)},        (F.i) 

where  the  various  functions  are  defined  as  follows 

M(,^)  =  61n^  +  lln(  -  iz(i)  -  Z(|)  -  l„c^  +  ?(f  -  ^)  , 

..2  1^2  'Ir^d-  C  1 


+  4c2Z(^)  +  ^^^^^  +  12Anc2  -  ^{1  +  2c2) 

^  s'^  2 

-3^[2Z(^)  +  4ln:^-5], 


2 


/_l(e,  A^)  =  61n^[l  +  2c4  -  4-^]  -  6Z(i)  -  12c4z(^)  -  12c'*lnc2 

.  '\/4  A  j2  Tyr2 

.Sa.2.V24;||,„|-2.Z,|^„, 


with 


^  r2^tan-l(lA4)  ,  (z  >  \) 

l.41n[(l  +  A)/(l-A)].     {z<\)    ,  (F.3) 

A=  11  -4z|2  . 
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APPENDIX  G 
'METAPLECTONS' 


G.l  Introduction 


Theories  of  objects  with  arbitrary  spin  and  statistics  known  as  anyons  are 
possible  in  (2  +  1  )-dimensions  [83].  Such  theories  have  been  apphed  suc- 
cessfully in  condensed  matter  physics.  For  example,  they  have  been  used  in 
the  description  of  quasiparticles  in  the  fractional  quantum  Hall  effect.  Some 
of  these  theories  employ  a  topological  quantity  known  as  the  Chern- Simons 
term.  This  term  is  sometimes  generated  radiatively  through  the  use  of  heavy 
fermions  rather  than  put  in  by  hand.  This  work  will  consider  such  theories. 
The  statistics  of  the  anyons  in  these  theories  will  depend  on  the  coefficient 
of  the  Chern-Simons  term  which  subsequently  will  be  seen  to  depend  on  the 
number  of  fermions  and  on  their  respective  charges.  If  all  the  fermions  carry 
the  same  fundamental  charge,  then  by  varying  the  number  of  fermions  a  cer- 
tain "spectrum"  of  fractional  statistics  is  possible.  This  spectrum  is,  however, 
limited  in  its  extent,  and  the  question  may  be  raised  as  to  whether  it  can  be 
broadened  by  effecting  changes  or  extensions  to  the  fermion  sector.  This  work 
addresses  this  question.  The  answer  is  yes,  and  the  change  or  extension  to  the 
model  considered  here  will  involve  fields  in  a  spinor  or  metaplectic  represen- 
tation of  the  orthosymplectic  Lie  supergroup.  The  representations  are  infinite 
dimensional  but  have  graded  dimension  that  is  finite  and  can  be  fractional. 
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In  Section  G.2,  a  brief  introduction  to  anyons  is  given  as  well  as  an  account 
of  the  role  of  the  Chern- Simons  term  in  such  theories.  The  specific  model  be- 
ing considered  is  presented,  and  the  extension  of  the  fermion  sector  is  moti- 
vated. Section  G.3  contains  a  self-contained  exposition  of  the  orthosymplectic 
Lie  supergroup  and  its  spinor  (or  singleton)  representations.  The  oscillator 
construction  is  discussed  and  used  to  arrive  at  the  metaplectic  representations. 
Section  G.4  applies  the  results  of  Section  G.3  in  extending  the  model  of  Section 
G.2.  The  "new"  theory  has  a  broader  spectrum.  The  statistics  are  tuned  in 
the  fermion  sector  now  extended  by  the  "metaplectons" .  A  short  summary  of 
the  results  of  the  appendix  and  some  remarks  make  up  the  contents  of  Section 
G.5. 


no 

G.2  Anyons 

In  (2  +  l)-dimensions,  the  faxniliar  Bose  and  Fermi  statistics  are  special 
cases  of  more  general  fractional  statistics.  Under  the  interchange  of  identical 
particles,  the  wave  function  changes  by  an  arbitrary  phase  rather  than  just 
being  symmetric  or  antisymmetric  as  in  3  +  1  dimensions.  An  anyon  is  any 
(2  +  1)- dimensional  object  with  fractional  statistics.  Such  an  object  can  be 
viewed  as  a  composite  consisting  of  a  "bare"  Bose  or  Fermi  point  charge  bound 
to  an  infinitely  thin  flux  tube.  Interchanging  two  of  these  identical  composites 
yields  any  phase  depending  on  how  the  fluxes  are  tuned.  The  phase  is,  in 
fact,  an  Aharanov-Bohm  phase  resulting  from  the  long  range  interaction  of 
one  composite's  charge  with  the  other  composite's  flux  tube.  A  paxticle  of 
charge,  Q,  encircling  an  infinitely  thin  flux  ($)  tube  picks  up  an  Aharanov- 
Bohm  phase,  e^'^-^B,  where  ^45  =  Q^    [84]. 

A  way  of  implementing  fractional  statistics  in  certain  (2  +  1)  models  is 
through  the  use  of  an  Abelian  Chern-Simons  term.  A  Lagrangian,  £,  with 
a  conserved  current,  j^,  can  have  additional  gauge  invariant  interactions  ap- 
pended to  it    [85] 

A£  =  -j^A^,  +  '^ae^'^'PA^duAp  ,  (G.2.1) 

where  A^  is  the  so-called  statistics  gauge  field  and  has  dimension  of  [ma55]2, 
and  a  is  the  Chern-Simons  coupling  with  dimension  of  [mass].  The  equation 
of  motion  of  this  gauge  field  gives, 

j^^  =  aet'^PduAp  .  (G.2.2) 

Integrating  the  zeroth  component  of  this  equation  over  all  space  yields 

I  d^xj"  =  a  I  d'^xB  ,     where    B  =  e'-''a.4,-  , 
J  J  ^  (G.2.3) 


Ill 

Hence  charged  particles  carry  statistical  flux  equal  to  Q/a.     This  result  is 
independent  of  the  insertion  of  a  kinetic  term  for  the  statistics  gauge  field 

^kin  =  -\F^"Ff,^  ,  (G.2.4) 

where 

Ff""  =  a^A'^  -  d^'A^'  .  (G.2.5) 

As  seen  in  the  following  argument,  if  a  kinetic  term  were  present,  the  equation 
of  motion  would  be  altered  to 

j^  =  aef^'Pd^Ap  +  duF^"'  .  (G.2.6) 

The  "statistics"'  photon  in  the  quantized  theory  will  have  a  mass,  a,  due  to 
the  Chern-Simons  term;  a  will  act  as  a  long  range  cut-off  for  the  electric 
and  magnetic  fields  [86].  Therefore  the  integral  over  all  space  of  the  time 
component  of  Eq.  (G.2.6)  will  still  give  Eq.  (G.2.3)  since  f  d'^xV  ■  'E  vanishes. 
In  a  Lagrangian  with  Eq.  (G.2.1)  added  to  it,  an  effective  current,  j*",,, 
can  be  defined  as  the  coefficient  of  d^x  under  the  gauge  transformation  A^  -^ 
A^^-d^X    [87] 

^eff  =  ^'  -  l^^'^'d^Ap  .  (G.2.7) 

Making  use  of  the  equation  of  motion,  Eq.  (G.2.2),  results  in 

Jeff  =  V^  ■  ^^-^-S) 

Accordingly,  this  leads  to  an  effective  charge  in  terms  of  the  physical  charge, 

Q, 

Qeff  =  ^Q  ■  (G.2.9) 

In  this  way  an  anyon  arising  from  a  Chern-Simons  term  can  be  considered  as 
point-particle-flux-tube  composite  with  statistical  charge,  Qeff,  and  flux,  $. 
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If  two  such  objects  are  adiabatically  interchanged  by  moving  each  one  half  way 
around  the  other,  the  Aharanov-Bohm  interchange  phase  induced  in  the  wave 

function  is 

^AB  =  oQeff^  +  oQeff^ 

-  -  (G.2.10) 

=  Qeff^- 
Inserting  the  results  of  Eqs.  (G.2.3)  and  (G.2.9)  above  yields 

Bab  =  ir  ■  (^-^-ii) 

The  full  quantum  statistics  phase  is  arrived  at  by  also  including  the  bare  statis- 
tics phase  {B\,art) 

Q^have+HB-  (G.2.12) 

Consider  the  following  functional  for  the  statistics  gauge  field,  A^J,,  in  (2  + 
l)-dimensions  involving  charge,  e  (which  has  dimension  of  [mas^js),  massive 
fermions  (which  have  dimension  of  [ma^^]) 

A[A]  =    /  Dil)Dli}  exp{i   /  (fix^^-f^^id^j,  -  eA^j,  -  m)V'} 

=  det(i:^-m)  , 
where  Tp  =  J^idfj,  +  ieAfx).  Through  vacuum  polarization  a  Chern-Simons  term 

is  induced  to  lowest  order  in  the  momentum    [88] 

lnA[A]  =  i^^  fd^xe^^PAf^d^Ap  +  I[A]  .  (G.2.14) 

Stt  |m|  y 

The  higher  order  terms  in  I[A]  are  proportional  to  powers  of  e  and  1/m,  and 

these  are  assumed  to  satisfy  the  condition 

m>e2.  (G.2.15) 

Comparing  expression  (G.2.14)  with  the  second  term  on  the  right  hand  side  of 
Eq.  (G.2.1)  and  assuming  the  fermion  masses  are  positive  yields  the  identifica- 
tion 

e2 

a=  — .  (G.2.16) 

Air 
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Note  that  the  induced  (Chern-Simons)  term  violates  parity,  a  fact  that  is  not 
surprising  since  a  fermion  mass  term  is  parity  violating  in  (2  +  l)-dimensions. 
Under  parity,  defined  by 

{t,x,y)  -^  {t,-x,y)  , 

(the  usual  definition  of  parity,  (x,y)  ->  (-x,  -y),  is  a  rotation  is  this  case)  the 
fermion  mass  term  transforms 

mtptp  — *•  —miptp  .  (G.2.17) 

Suppose  there  had  been  A''  species  of  fermions  with  charges  e^,  k  =  1,...,N, 
then 

f  -  /•  ^  _ 

ifc-i 

^  (G.2.18) 

=  H  detii-pk  -  m)  , 
k=l 

where  pk  =  7/^(5^  +  iCkA^).  The  induced  Chern-Simons  term  then  yields  the 

following  identification  for  the  coupling 

1     ^ 
''  =  ^T.4-  (G.2.19) 

k=i 

Some  models  exhibiting  fractional  statistics  and  employing  a  Chern-Simons 
term  sometimes  replace  it  with  a  Lagrangian, 

^N  =  J2  ^-^-(^^  ~  e4-  m)V't  ,  (G.2.20) 

k=l 

thereby  generating  the  Chern-Simons  term  radiatively  as  previously  discussed. 

This  will  yield  a  Chern-Simons  coefficient  given  by  Eq.  (G.2.19)  with  e^  =  e 

for  all  k  =  l,...,iV 

«  =  N—  .  (G.2.21) 

47r  ^  ' 
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Consider  the  following  Lagrangian    [89] 

£  =  |(a^  +  iqeA^)<p\'^  -  A/2|<;5|2  -  ^F'^^'F^^  +  Cn  •  (G.2.22) 

The  bare  boson  (Oi^re  —  0)  field,  <^,  carries  charge  Q  =  ge  (g  is  an  integer) 
and  has  mass,  M.  Interchanging  two  of  these  quanta  will  give  Aharanov-Bohm 
phase  (G.2.11).  Hence  their  quantum  statistics  are  given  by 

^=-^.  (G.2.23) 

If,  for  example,  q  =  I  and  N  =  2,  then  the  statistics  are  those  of  fermions 
[0  =  tt).  More  exotic  phases  are  possible  with  more  fermions  (i.e..  N  >  3). 
Indeed,  for  q  —  I  and  N  =  4,  the  statistics  are  those  of  semions  (6  =  7r/2). 

It  is  worth  noting  that  a  composite  of  N  ajiyons  has  an  interchange  phase 
i^ABi^))  that  is  N  times  the  interchange  phase  of  one  anyon  [90].  Therefore 
in  this  case  the  full  quantum  statistics  phase  is 

(G.2.24) 

=  ^harei^)  mod  2tt  . 
Hence  such  an  object  (an  N-anyon)  exhibits  no  fractional  statistics  but  rather 
exhibits  its  bare  statistics. 

Evidently  the  Lagremgian  (G.2.22)  which  has  a  global  U{N)  symmetry 
cannot  produce  arbitrary  fractional  statistics  for  a  given  value  of  q,  rather  the 
statistics  phases  are  constrained  by  the  integer  ^V  through  Eq.  (G.2.23).  One 
way  of  extending  the  range  of  0  is  by  having  non-integer  N .  This  may  seem 
impossible  if  N  is  thought  of  as  the  number  of  fermions  but  not  so  impossible  if 
it  is  identified  with  the  power  of  the  determinant  (G.2.18).  The  purpose  of  this 
work  is  to  point  out  that  non-integer  powers  of  this  determinant  are  attainable 
by  considering  more  exotic  internal  symmetry  groups.     These  produce  new 
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types  of  fractional  statistics,  in  the  sense  discussed  above,  hitherto  unknown  in 
the  Hterature.  This  is  introduced  through  the  use  of  supergroup  representations 
with  fractional  dimensions. 

G.3  The  Metaplectic  Representations  of  0Sp{r/2m,R) 
The  use  of  supergroup  representations  as  a  device  to  yield  fractional  powers 
of  determinants  is  not  unknown  in  the  literature  [91].  The  covariant  quan- 
tization of  the  Green-Schwarz  superstring  appears  to  require  the  introduction 
of  infinite  sets  of  extra  fields,  of  ghosts,  and  of  ghosts  for  ghosts.  It  has  been 
suggested  that  the  superstring  action  can  be  written  in  an  OSp  invariant  form 
for  certain  choices  of  gauge  with  the  "infinite  towers"  of  ghosts  transforming  as 
spinor  (or  metaplectic)  representations  of  the  0Sp{r/2m,R)  Lie  supergroups 
[92].  The  haxmonic  oscillator  method  for  constructing  unitary  lowest  weight 
representations  (i.e.,  representations  in  which  one  of  the  generators  is  bounded 
from  below)  of  the  noncompact  supergroup,  OSp{2n/2m,R),  has  been  pre- 
sented in  the  literature  [93].  In  the  following  the  reader  is  assumed  to  have 
some  familiarity  with  Lie  supergroups  and  superalgebras  (for  a  review  see  ref- 
erence   [94]). 

The  orthosymplectic  Lie  supergroup,  OSp(2n/2m.  i?),  has  even  subgroup, 
S0(2n)  X  Sp{2m,R).  The  number  of  even  generators  is  ^(2n)(2n  —  1)  + 
5(2m)(2m  -f  1).  The  number  of  odd  generators  is  (2n)(2m).  OSp{2n/2m,R) 
has  li(m/n)  as  its  maximal  rank,  compact  subsupergroup.  With  respect  to 
this  subsupergroup,  OSp{2n/2m,  R)  admits  a  Jordan  decomposition  or  three- 
grading 

L  =  L.i®Lo^L+,  .  (G.3.1) 
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L  and  Lq  denote  the  Lie  superaigebras  of  the  supergroups  OSp{2n/2m,  R)  and 
U{m/n),  respectively.  U{m/n)  has  even  subgroup  U{n)  x  U{m)  with  respect 
to  which  S0{2n)  x  Sp{2m,R)  has  a  three-grading.  The  Lie  superalgebra  of 
OSp(2n/2m,R)  is  reahzed  as  bilineaxs  of  a  set  of  /  =  2p  +  £  superoscillators 
(and  their  hermitian  conjugates) 

(^(s)  =  ei(^)  ,    V^{s)  =  r,\is)  ,    C"^  =  Ci  (G.3.2) 

i  =  l,...,m;    /i  =  l,...,n;    5  =  1,..., p. 

Depending  on  whether  there  is  an  even  or  odd  number  of  superoscillators,  e 
will  equal  0  or  1.  The  metaplectic  (or  singleton)  representations  correspond 
to  /  =  1  {p  —  0,e  =  1).  The  first  m  components  of  the  superoscillators 
are  bosonic  while  the  last  n  components  are  fermionic.  The  superoscillators 
^A^  VA^  Ca  (s  1  ^  1  C  )  transform  in  the  covariant  ( contravariant )  fundamental 
representation  of  U{m/n)  and  satisfy  the  supercommutation  relations 

[a,c^}  =  <^i.  (G.3.3) 

All  other  supercommutators  vanish.  The  graded  commutator  is  defined  as 

lU^^}  =  U^^  -  {-l)'^^^'^^k^U  ,  (G.3.4) 

where  g{A)  is  0  or  1  depending  on  the  grade  of  the  index  .4  =  {i,iJ.).  L-i 
consists  of  di-annihilation  operators  of  the  form 

(a  •  nB  +  nA-(B+^  (aCb  ■  (G.3.5) 
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Correspondingly,  L+i  consists  of  di-creation  operators 

The  bilinears  of  Lo  are  of  the  form 

f^  •  (b  +  {-i)'^^^'^^Ub  •  ^-^  +  J  [C^C5  +  (-D^^^^^^^VbC-^I  ■      (G.3.7) 

The  vector  notation  employed  above  is  defined  as  follows 

P 
L-^B^J2^A^^^'Jb{^)  ■  (G.3.8) 

5=1 

The  unitary  lowest  weight  representations  of  OSp  are  constructed  by  consid- 
ering a  set  of  states,  1^2),  in  the  (super)  Fock  space  of  superoscillators.  These 
states  transform  irreducibly  under  Lo  (the  Lie  algebra  of  the  maximal  compact 
subsupergroup)  and  are  annihilated  by  the  operators  belonging  to  L_i.  Op- 
erating successively  on  \Q.)  with  the  elements  of  X+i,  a  unitary  lowest  weight 
representation  is  generated 

{|fi)  ®  L+ilQ)  e  (L+O^lfi)  ©■••}.  (G.3.9) 

For  /  =  1,  \Q.)  =  {lO),^  |0)}  are  the  only  two  lowest  weight  states,  where  |0) 
is  the  Fock  space  vacuum  annihilated  by  all  annihilation  operators. 

Application  of  the  oscillator  construction  to  S0(2n )  or  Sp(2m,  R)  is  straight- 
forward. In  the  case  of  S0(2n),  the  set  of  /  =  2p+£  fermionic  oscillators  would 
consist  of 

a/i(-s)  ,    l^fi{s)  ,    7/i 

a^(s)  =  al(s),    ;3^'is)  =  0l(s)  ,    1^  =  ll  (G.3.10) 

f^i  =  1 n  ;    5  =  1,  ...,p  . 
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Clearly  the  form  of  the  unitary  irreducible  lowest  weight  representation  (Eq.  (G.3.9)) 
indicates  that  it  is  finite  dimensional  when  the  oscillators  axe  fermionic,  as  is 
the  case  for  S0{2n).  For  Sp{2m,R),  the  oscillators  are  bosonic,  and  the  rep- 
resentations are  infinite  dimensional.  The  unitary  irreducible  representations 
can  be  put  in  a  form  that  displays  the  even  subgroup  (S0(2n)  x  Sp{2m,R)) 
structure.  The  two  irreducible  spinor  representations  of  OSp{2n/2m,  R)  can 
then  be  written 


*+  = 


/(2) 


,  *-  = 


/^(0)\ 

V  '■■  J 


(G.3.11) 


The  V''s  are  SO{2n)  chiral  spinors  with  the  Grassmann  parity  alternating  down 
each  "tower"  (the  Grassmann  parity  of  the  first  in  the  sequence  is  arbitrary). 
The  superscript  indicates  their  transformation  properties  under  U{m),  the 
maximal  compact  subgroup  of  Sp{2m,R).  For  example,  rp_^  =  0^"'*  is  a 
(-f  )-chirality  S0{2n)  spinor,  Grassmann  odd  field  transforming  as  a  ^  rank 
symmetric  tensor  of  SU{m).  It  is  possible  to  combine  the  two  irreducible 
representations  to  form  a  "Dirac"'  OSp  spinor 


*  = 


/^(0)\ 
^(1) 
^(2) 


V' 


(") 


(G.3.12) 

V  ■■  ) 

The  two  irreducible  ((  +  )  and  (  — )  chirality)  spinors  (in  the  basis  that  they  have 
been  presented  above)  can  be  recovered  by  applying  the  following  projection 
operators 

P±  =  1®P±,  (G.3.13) 


where  1  is  the  infinite  unit  matrix  and  P±  are  the  S0(2n)  projectors 


1 


p±  =  -(i±r),   r 


1       0 
0     -1 


(G.3.14) 
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and  r  has  been  displayed  in  2  x  2  block  form.  50(2n  + 1)  has  no  such  projectors 
and  correspondingly  0Sp{2n  +  l/2m,  i?)  has  only  one  irreducible  singleton 
representation.  Indeed,  instead  of  a  three-grading,  0Sp{2n  +  l/2m,R)  admits 
a  five-grading  (Kantor  decomposition)  with  respect  to  its  maximal  compact 
subsupergroup    [95] 

I  =  L_i  ®L_i  ®io®i+^  ®i+i  •  (G.3.15) 

The  invariant  inner  product  of  complex  OSp  spinors  is  given  by 

l^A  =  J2^''"'^''"''''^idr--^hh  '  (G.3.16) 

where  ^  and  A  are  two  OSp  spinors,  ^  =  '^'Tt  is  the  Pauli  adjoint  and  Ft 
represents  the  product  of  the  time-like  gamma  matrices. 

The  graded  dimension  of  a  representation  is  defined  as  the  number  of  Bose 
states  minus  the  number  of  Fermi  states  in  the  vector  space.  It  can  be  computed 
by  evaluating  the  character  for  the  identity  matrix  [96].  Consider  an  irreducible 
spinor  representation  of  0Sp{r/2m,R).  If  r  =  2n,  then  the  dimension  of 
an  50(2n)  chiral  spinor  is  dspinoriSO{2n))  =  2""^.  If  r  =  2n  +  I,  then 
dspinor{SO(2n  -f-  1))  =  2".  The  dimension  of  a  k  rank  symmetric  tensor  of 
SU{m)  is  given  by  the  following  binomial  coefficient 

dsym.k-tensoA^Uim))  =  r" ^  _  _^      j.  (G.3.17) 

If  each  Grassmann  odd  field  counts  as  —1  and  each  Grassmann  even  field  counts 
as  -f-l,  then  the  following  regulated  sum  of  the  SU{m)  symmetric  tensors  in  ^ 
yields 

^^^  V     m-1     ;        .-.if^^  \     m-l     J 

=  lim(l  +x)~'" 

X 1 

.-y—m 

(G.3.18) 


120 
The  graded  dimension  of  an  0Sp{r/'2m,  R)  spinor  is  therefore  given  by    [97] 

dmetaiOSp{r/2m.R))  =  2-"^J,p,w(50(r))  .  (G.3.19) 

G.4  New  Representations  of  Anyons 
The  supergroup  representations  of  the  last  section  will  now  be  used  to  show 
how  to  extend  the  statistics  spectrum  discussed  previously.    As  an  example, 
consider  the  following  functional  of  A^ 

/oo  .00 

Yl  D^l^kDi^k  exp{i   /  d^x  Yl  ^k(i^  -  e4-m)^k} 
k=i  •'  fc=l  (G.4.1) 

=  [det(z:^-m)](^-^+^-l+-)  , 
where  the  V'/fc's  are  (2  +  1)  space-time  spinors  having  alternating  Grassmann 
parity  with  the  odd  subscript  ones  odd  and  the  even  subscript  ones  even.  This 
sequence  of  fields  {V'/fc}^  ^  =  1^2, ...,  can  be  fitted  into  a  metaplectic  represen- 
tation of  05p(2/2,  R) 


*  =  U..        •  (G.4.2) 


The  V'jt's  are  one  component  chiral  50(2)  spinors  and  are  scaiars  under  5p(2,  R) 
(a  mass  term  is  possible  for  chiral  50(2)  spinors).   The  Lagrangian  for  these 

fields  can  be  rewritten  using  the  boldface  notation  introduced  in  Section  G.3 

00 

Cmeta  =  ^{ip-  m)^  =  J^  ^ki^f-  "^)^k  ■  (G.4.3) 

k=l 

The  alternating  sum  in  the  power  of  the  determinant  of  Eq.  ( G.4.1 )  can  be  com- 
puted in  a  natural  way  using  the  group  invariant  Introduced  in  Eq.  (G.3.19), 

dmetaiOSp(2/2,R))^l.  (G.4.4) 
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Therefore, 

^meta[A]  =  [det(2:^-  m)]2   .  (G.4.5) 

The  induced  Chern-Simons  term  will  then  have  coupHng, 

The  square  root  of  the  Dirac  operator  enters  the  analysis  of  the  SU{'2) 
global  anomaly  in  four  dimensions  [98].  There,  integrating  over  an  SU{2) 
Weyl  doublet,  rp,  gives 

I DxbD^  exp{i  I  d^x^fip}  =  [detz:^]2  .  (G.4.7) 

Such  a  representation  for  the  Weyl  fermions  is  possible  because  SU{2)  has  no 
local  anomaly 

Tr({T'',T''},T')  =  0  ,  (0.4.8) 

where  the  T's  are  SU[2)  generators.  In  the  present  case  no  such  condition 
need  be  satisfied  by  the  gauge  group.  The  nontriviality  of  the  homotopy  group, 
Il4{SU(2)),  is  crucial  in  the  analysis  of  the  global  anomaly.  In  three  dimensions 
the  relevant  group  is  Eg.  However,  Il3{U{l))  =  0.  and  hence  no  global  anomaly 
plagues  the  metaplecton  theory  being  considered  here.  Clearly  metaplecton 
theories  involving  non-Abelian  gauge  groups  with  nontrivial  homotopy  groups 
may  suffer  from  such  problems. 

Using  the  result  of  Eq.  (G.4.6),  the  Lagrangian, 

C  =  \{d^  +  igeA^)<^|2  -  Af2|<^|2  _  ^Ff'^F^,^  +  Cmeta  ,  (G.4.9) 

will  lead  to  a  statistics  phase, 

d  =  iwq'^,  (G.4.10) 
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for  the  charge  qe  excitations;  in  other  words,  they  are  bosons.  It  appears  that 
the  OSp(2/2,  R)  global  symmetry  has  alone,  at  least,  not  produced  any  differ- 
ent statistics  (in  the  sense  of  Section  G.2).  However,  consider  the  Lagrangian, 

£'  =  \{^^,  +  iqeAi,)<t>\^  -  M2|,^|2  -  ^F^^F^,^  +  Cmeta  +  C^'  ,  (G.4.11) 

where 

N'  _ 

£^,  =  Y,  n{^'^  -  e-A-  m)0jt  •  (G.4.12) 

This  model  has  a  U(N')x  05p( 2/2,  R)  global  symmetry  and  will  lead  to  Chern- 
Simons  coupling, 


2N'  +  1,  e^ 
a  = 

and  to  statistics  phase, 


a=^^— (-),  (G.4.13) 


^  =  :T^  •  (G.4.14) 

2N'  +  1  ^  ^ 

This  phase  (taking  N'  >  1}  cannot  be  achieved  in  model  (G.2. 22)  for  any  value 
of  N  {N  an  integer).  So  the  goal  of  attaining  a  wider  range  of  fractional  statis- 
tics has  been  reached  through  the  addition  of  Cmeta-  More  phases  outside  the 
range  of  (G.2. 22)  are  possible  by  generalizing  to  OSp(2/2m,  R),  with  the  fields 
in  Cmeta  transforming  as  the  metaplectic  representation  of  this  supergroup. 
The  statistics  phase  of  a  U{N')  x  OSp{2/2m,  R)  model  would  be 

e=——^ ,  (G.4.15) 

N'  +  dmetairn) 

where  (using  Eq.  (G.3.19)) 

dmetaim)  =  dmetaiOSp{2/2m,  R))  =  2-"*  .  (G.4.16) 

It  should  be  emphasized  that  the  statement  that  Eq.  (G.4.14),  associated 
with  Lagrangian  (G.4.11),  gives  phases  outside  the  range  of  (G.2. 22),  is  a  valid 
one  for  fixed  q  in  (G.2. 22)  and  (G.4.11).  Indeed  statistics  (G.4.15)  with  given 
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values  of  q,  N' ,  and  m  are  attainable  using  (G.2.22)  with  no  OSp  symmetry. 
However,  the  (f)  field  would  then  have  statistics  charge,  q' ,  such  that 

q'  =  \'^ll^        "^^"^"'  (G.4.17) 

I  225  ,     m  odd    ,  ^ 

and  there  would  have  to  be  more  than  N'  fermions 

iV-/2'"iV'  +  l,  meven    ;  (CA^R^ 

^^  -  1  2(2'".V' +  1)  ,     modd    .  (G.4.18) 


G.5  Conclusion 

This  appendix  has  shown  how  to  generate  a  new  class  of  models  of  anyons  in 
(2  +  1  )-dimensions  employing  a  radiatively  generated  Chern-Simons  term  that 
have  their  "spectrum"  of  fractional  statistics  broadened  without  varying  any  of 
their  parameters  [11].  The  fields  introduced  to  induce  the  Chern-Simons  term 
are  supplemented  with  fields  forming  a  metaplectic  or  spinor  representation 
of  the  Lie  supergroup,  05p(r/2m,  R).  The  representations  are  infinite  dimen- 
sional, since  the  even  Lie  subgroup  of  0Sp{r/2m,  R)  is  SO{r)  x  5p(2m,  R),  and 
Sp{2m,  R)  has  infinite  dimensional  spinor  representations.  These  metaplectic 
representations  have  the  special  property  that  they  yield  fractional  powers  of 
determinants  when  the  fields  are  integrated  out.  In  the  models  mentioned 
above  they  lead  to  statistics  outside  the  range  achievable  with  any  number  of 
ordinary  fermions. 

There  are  many  ways  to  generate  exotic  statistics  in  (2  -|-  l)-dimensions. 
This  work  has  presented  a  new  method.  In  general  the  distinguishing  feature 
between  these  methods  is  the  fermion  representation.  As  it  stands,  the  new 
fermion  representations  introduced  are  indistinguishable  from  any  other  since 
the  fermions  are  integrated  out.  In  this  regard  the  new  representations  are 
nothing  more  than  a  device.    To  distinguish  two  such  equivalent  theories,  it 
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is  necessary  to  either  include  more  degrees  of  freedom  such  as  scalars  trans- 
forming nontrivially  under  the  symmetry  group  or  to  gauge  some  parts  of  the 
group.  Whether  or  not  this  will  lead  to  any  interesting  results  is  not  known. 
However,  as  it  stands  the  model  "forgets"  about  the  metaplectons  after  they 
are  integrated  out.  Gauging  is  a  means  of  effectively  leaving  a  "memory"  of 
the  OSp  group  after  the  metaplectons  are  integrated  out. 

This  appendix  has  not  dealt  with  the  dynamics  of  anyons  but  rather  with 
an  alternate  means  of  inducing  fractional  statistics.  It  is  noted  that  questions 
of  this  sort  have  been  addressed  elsewhere,  and  a  relativistic  wave  equation  for 
anyons  has  been  formulated  derivable  from  an  appropriate  action    [99]. 
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